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Abstract 

A method is proposed to calculate the confidence level for agreement of data 
with the Standard Model (SM) by combining information from direct and indirect 
Higgs Boson searches. Good agreement with the SM is found for niH — 120 GeV 
using the observables most sensitive to tjih- Ai and mw- In particular, quantum 
corrections, as predicted by the SM, are observed with a statistical significance 
of forty-four standard deviations. However, apparent deviations from the SM of 
3.7cj and 2.8cr are found for the Zz^z7 and right-handed Zbb couplings respectively. 
The maximum confidence level for agreement with the SM of the entire data set 
considered is ~ 0.006 for mn — 180 GeV. The reason why confidence levels about 
an order of magnitude higher than this have been claimed for global fits to similar 
data sets is explained. 



1 Introduction 



It is now almost three decades since the the first accelerator [1, 2] and physicss [1, 3] 
studies, that eventually lead to the construction and operation of the LEP e+e~ collider 
at CERN were performed. Now, some roughly twenty man-millenia of work by physicists 
and engineers later, the almost final results of the LEP and concurrent SLC (Stanford 
Linear Collider) experimental programs are available [4, 5]. By general consensus, the 
most scientifically important of these results concern high precision tests of the Standard 
Electroweak Model (SM) [6, 7, 8]. During the same period, important contributions to this 
subject (discovery of, and measurement of the mass, nit, of the top quark, measurement 
of the mass, mw, of the W boson and the NuTeV neutrino-quark scattering results) were 
also made at the FERMILAB laboratory. 

The LEP program consisted of two stages. In the first, 'Z-pole', running a total of 
~ 1.7 X 10^ Z decays into pairs of all the fundamental fermions (the matter fields of the 
SM), except the top quark, were collected by the four LEP detectors, ALEPH, DELPHI, 
L3 and OPAL. At SLC, the same processes were studied with lower statistics but with 
precisely controlled electron beam polarisation, allowing experiments of very high sesitiv- 
ity to be performed, so producing results of comparable statistical accuracy to those of 
LEP. The final outcome of all this work may be summarised, in what concerns the SM, in 
only seven numbers, which, for definiteness can be taken to be the right- and left-handed 
couplings to the Z of the charged lepton, c-quark and b-quark pairs and the (left-handed) 
coupling of the Z to neutrino pairs. The busy reader, who would like to go straight to the 
final conclusions, is invited to look directly at Tables 19 and 20 below where the measured 
values of these coupling constants, together with the corresponding SM predictions are 
shown. The second, 'high energy' phase of LEP operation at collision energies above the 
threshold for W pair production provided essentially one additional high precision SM 
parameter, mw This measurement, in combination with the FERMILAB measurement 
of comparable accuracy of the same quantity, is also shown in Tables 19 and 20. Other 
important measurements performed during the second LEP phase were of the triple boson 
WW7 and WWZ couplings, but these have less impact than the fermion couplings as a 
test of the core physics underlying the SM: a renormalisable quantum field theory incor- 
porating local gauge invariance that is spontaneously broken by the Higgs mechanism^. 
They are therefore not discussed further in this paper. In the conventionally used on- 
shell renormalisation scheme [9] where niw is traded as an input parameter for the much 
more precisely known Fermi constant, G^, derived from the measured muon lifetime, 
the predictions of the SM depend (apart from other and better known parameters) on 
three relatively poorly known ones: rrit, the electromagnetic couplant at the Z mass scale, 
a{mz), and the mass of the Higgs boson, niH- Because of quantum loop corrections, the 
right- and left-handed coupfings, and , of the Z to charged lepton pairs, as well as 
mw, are strongly sensitive to the values of mt and mn- Actually, it is only the ratio 
Qi'/gi and mw which are strongly sensitive to mn, so that the measurements of these 

-'^In fact the now experimentally well verified [4] SM predictions for the WW7 and WWZ couplings 
were aleady contained in Glashow's original electroweak paper [6] where they follow at tree level from 
global SU(2)i:, invariance and quantum mechanical mixing of the W° and B*^ fields. Their values then 
shed little light on the correctness (or otherwise) of either the rcnormalisabilty of the theory or the Higgs 
mechanism. A genuine test of local gauge symmetry would be provided by measurements of the strength 
of quadrilinear boson couplings. So far this has not been done. 
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two quantites provide the most stringent limits, from quantum corrections, on the value 
of niH- 

In the second phase of LEP, an unsuccessful search was performed for directly produced 
Higgs bosons [10], resulting in the 95% confidence level (CL) lower limit uih > 114.4 GeV. 
The principle aim of the present paper is to combine, in a transparent way, this direct 
limit with the indirect information derived from 'g^/gi (or equivalently Ai, see below) 
and rriw, to derive combined curves of CL, the confidence level that the mi:/-sensitive 
data agrees with the SM, as a function of rriH- The reader might then reasonably hope 
that the result of the paper would be a single curve of CL versus ttih- In fact she (or 
he) will find eight figures where CL is plotted versus ttih containing in total 18 different 
curves. The reason for this complication is that the data, though perfectly consistent 
exprimentally^, is not consistent with the SM (see Tables 19 and 20) and depending on 
the assumptions made (SM correct, model-independent analysis, certain data included 
or excluded) different results are found for the CL curves. I have included a number of 
different possibilities to demonstrate the inconsistency of the data with SM predictions. 
The reader may then choose the curve for which the assumptions match best her (or his) 
own favourite ones. 

The structure of the paper is as follows: In the next Section some general remarks 
concerning the different functions of the science of Statistics in data analysis are made. In 
particular it is pointed out that the choices made until now for the estimator in global 
electroweak analyses give an over-optimistic estimate of the level of agreement of the data 
with the SM. In Section 3, the heavy quark asymmetry measurements that are prima 
facie inconsistent with charged lepton asymmetries, when both are interpreted within 
the SM, are discussed. In particular, the internal consistency of the data and systematic 
error estimates are examined in some detail. The sensitivities of different electroweak 
observables to rrit and rriH are discussed in Section 4. This is the only place in the paper 
where fit results are shown and discussed. It is demonstrated that the sensitivity to uih 
comes essentially from only the observables Ai and mw- Section 5 describes the algorithm 
used for combining the direct and indirect measurements of ttih- Results for CL derived 
from Ai and mw, assuming the correctness of the SM, but selecting different data, are 
shown. Also shown in this Section is the sensitivity of CL to the values of the parameters 
rrit and a{mz)- In Section 6 the alternative interpretations of the result of the NuTeV 
experiment are explained and it is pointed out that the interpretation, as required in a 
model-independent analysis, as a measurement of the Zi/V coupling, instead of mw, is 
strongly favoured by arguments of statistical consistency. In Section 7 a complete set 
of model independent observables is extracted and compared with SM predictions. Con- 
straints are set on the coupling of non-b down-type quarks using the precisely measured 
observable Thad- Quantum corrections are extracted for different fermion flavours and 
compared with SM predictions. Finally, in Section 7, curves of CL versus mn derived 
from a estimator using all or selected subsets of the considered observables (including 
now mj:/-insensitive ones) are shown. In Section 8 values of CL obtained as described in 
Section 7, are compared with the confidence levels of previously published global fits to 
similar data. The confidence levels are seen to be very consistent when the purely sta- 
tistical dilution of the hypothesis testing power of the estimators of the global fits, as 

^That is, good agreement is found between different measurements of the same experimental observ- 
ables. For details see Reference [4] 
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discussed in Section 2, is taken into account. Section 8 also contains a critical discussion 
of niH limits determined from Ax^ plots. Section 9 contains a summary and conclusions, 
including the author's personal choice of three most pertinent CL versus rriH plots. The 
busy reader is encouraged to read this Section first to get a general view of the results and 
conclusions returning later (if still interested) to the earlier Sections for more information 
and supporting arguments. 

When the first version of the present paper was almost complete a new experimental 
value of rrit, 178 ± 4.3 GeV, was announced by the CDF and DO collaborations [11]. 
Since the change from the previous value of 174.3 ±5.1 GeV has a dramatic effect on the 
CL curves, especially for large values of ttih, all such curves shown, when the contrary 
is not explicitly stated, use the new value of rrit. However no global fits have yet been 
published by the EWWG and EWPDG using the new value. Therefore, in Section 4 
where comparisons with the results of EWWG global fits are made, the old value of rrit is 
used. The conclusions of Section 8, where the consistency of the confidence levels found 
in the present paper with those quoted for global fits is discussed, are unaffected by the 
change in the measured value of rrit. 

2 Statistics: Data Consistency versus Hypothesis Test- 
ing 

In the context of the analysis of experimental data. Statistics has three quite distinct 
roles to play. These are: 

(i) To judge whether different measurements of the same physical quantity are con- 
sistent with each other, and to derive an unbiased weighted average value of the 
quantity. 

(ii) To test the hypothesis that an ensemble of measurements of the same or different 
physical quantities are consistent with some theory. 

(iii) In the case of positive answers to the questions implicit in (i) and (ii), to determine 
numerical values of unknown or partially known parameters of a theory from the 
data. 

In previous and current analyses of precision electroweak data performed by the LEP 
and SLD electroweak working groups (EWWG) [4] and the standard model sub-group 
of the Particle Data Group (EWPDG) [5], only the functions (i) and (iii) above are 
systematically performed, with little, if any regard for (ii). In fact tests of data consistency 
(comparisons of different measurements of the same physical quantity) are performed by 
the EWWG using the estimator, with, in general, very satisfactory results [4]. In the 
global fit to all data, since the point (ii) is not addressed, all relevant data is used for 
parameter estimation in the global fit. In the case that all of the data is in agreement with 
the SM this procedure gives the best, unbiased, estimate of parameter values. However, 
if certain sub-sets of data do not agree with the SM, biased results may be obtained 
using this procedure. In particular, as will be discussed below, the fitted value of mn 
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obtained with the current data set is biased towards higher values by about 50 GeV 
by just such an effect. The EWPDG also do not investigate the level of agreement of 
different data sub-sets with the SM and related biases, being concerned only with the 
function (iii), parameter estimation on the assumption that all data is correctly described 
by the SM [12]. In this case different measuremnts of the same physical quantity are 
included as independent data in the fit without any prior consistency checks such as those 
prrformed by the EWWG. In the case that subsets of data do not agree with the SM, 
fitted prameters may then be biased in just the same way as in the EWWG global fits. 
It seems to the present writer that seeking the answer to the question posed in (ii) above 
should be the principle aim of experimental investigations of the SM, but, as a point of 
fact, this is an avowed goal of neither the EWWG nor the EWPDG. 

So what is the answer to the question implicit in (ii) above provided by the current 
electroweak data set? The nature of the problem is well illustrated by some fit results 
quoted in a paper devoted to a search for possible evidence of supersymmetry in precision 
electroweak data [13]. Fitting, as a preamble, the minimal electroweak standard model 
to only the sin^ ®e// values derived from either leptonic or hadronic asymmetry measure- 
ments, a per degree of freedom (x^/d.o.f.) of 18.4/4 was obtained corresponding to 
a CL of 0.001. A fit by the EWWG to the same data set, but using instead about 20 
observables^ reported a x^/d.o.f. of 26.0/15, with a CL of 0.04. An analysis of essentially 
the same year 2000 data set by EWPDG, but fitting the SM to more than 40 observables 
found, for a global fit, a xVd-o.f. of 42/37 with a CL of 0.27 [14]. The fitted value of 
m//, was very similar in these three different fits since, as discussed below, almost all the 
sensitivity to vhh is found in only two observables, sin^ ©e// ^"^^ 'mw- Thus, for essentially 
the same fitted value of uih, CLs differing by a factor of up to 270, according to the fit 
procedure used, were obtained. Which (if any) of the different CLs most truly refiects the 
agreement between the data and the SM prediction? The principal aims of the present 
paper are, firstly, to provide an answer to this question, and, secondly, to combine CLs 
derived from direct and indirect experimental limits on itlh so has to obtain an meaning- 
ful overall CL that reflects both the internal consistency of different observables and the 
global level of agreement with the SM. 

The explanation of the poor CL obtained in the fit to only the leptonic and hadronic 
sin^ ©i// values is now well known. As first pointed out in analyses of the 1996 data 
set [15, 16] The Z-boson b-quark couplings appear to be anomalous at about the three 
standard deviation level. These couplings are quite insensitive, in the SM, to uih and rrit, 
but, due to a correlation effect, when heavy quark forward/backward asymmetries are 
analysed, assuming the correctness of the SM, to extract a value of sin^ ®e//' ^^'^ latter is 
found to correspond to a much larger value of uih than that derived from purely leptonic 
measurements [17, 18]. This leads to barely compatible values of sin^ Ggjy from leptonic 
and hadronic (essentially b-quark) data and explains the poor CL of the fit to this data 
to obtain ttih and rrit mentioned above. 

More recently, much more precise experimental measurements of mw have become 
available. These are found to favour a value of uih even lower than that suggested by the 

^Many of these quantities are actually 'pseudo-observables', but for brevity the term 'observable' will 
be used throughout this paper for extracted physical quantities sensitive to parameters of the electroweak 
theory. 
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leptonic data, thus resulting in a large discrepancy between the uih value obtained by 
combining the leptonic data and mw and that derived from hadronic asymmetries. This 
problem has been has been recently stressed in the literature [19] and is now generally 
appreciated [20]. 

The reason for the factor ~ 300 difference in the CLs of different fits is easily under- 
stood. The point is that the fit to only the sin^ values was essentally performing 
the function (ii) above, i.e. hypothesis testing, whereas the EWWG and EWPDG fits 
were combining the functions (i) and (ii) with a large weighting factor in favour of (i). 
How this happens will now be explained. In addition to this effect, the hypothesis testing 
abihty of the fit is further blunted by the inclusion of observables in the fit, that have 
almost no sensitivity to ttih and m^, in both the EWWG and EWPDG analyses. 



Consider a number, A^, of independent measurements, Qi, of the same quantity, Q. The 
theoretical expectation for Q is Qrhy and the weighted everage value of the measurements 
is Q. With the assumption of uncorrelated experimental errors, three different Pearson 

estimators may be defined, as follows: 

2-^ ^2 



Xdata.WA 



Xdata,Thy 



{Qi ~ QrhyY 

1=1 



2 _ {Q — QrhyY 

XwA,Thy - -2 



(2.1) 

(2.2) 
(2.3) 



In Eqn(2.3), a is the weighted mean error on the quantity Q. Assuming uncorrelated, 
Gaussian distributed, errors it is given by the relation: 



1 



N 

E 



1 



^ ^i 



where CTj is the estimated RMS uncertainty on Qi. Noting the identity: 

Qi — Qrhy = {Qi — Q) + {Q — Qrhy) 

Eqn(2.2) may be written as: 



(2.4) 
(2.5) 



X.data,Thy 



N 

E 
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N 
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N 
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i=i n 



U=l 

Q 



2^i = l <^2 



=1 n 



_ 2 I 2 

— Xdata,WA + XwA,Thy 



(2.6) 



where, in the third line of Eqn(2.6) the definition of ct, (2.4) and Eqn(2.3) have been 
used, and in the fourth line the definition of Q: 
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So, in the simple case of uncorrelated Gaussian errors, the for consistency of the 
data with the theory is equal to the simple sum of the for consistency of the data 
with its weighted average plus the for consistency of the theory with the weighted 
average. Clearly xlata WA ^ measure only of the internal consistency of the data and so 
the corresponding CL provides an answer only to the question raised in point (i) above. 
"XswA.Thy gives, providing the CL for ')dd.ata,WA acceptable, an estimate of probability 
that the data is correctly described by the theory, and so provides the hypothesis test 
mentioned in (ii) above, as well as estimating the values of unknown parameters of the 
theory (for example, mn if Thy = SM) in accordance with point (iii) above. However, 
the value of Xdata,Thy^ statistical estimator universally used by both the EWWG and 
the EWPDG, reflects both the internal consistency of the data and the level of agreement 
of the data with the theory. If the number of data is very large, the relative contribution 
of XwA,Thy to Xdata,Thy bccomcs Very Small, since the former has only one degree of 
freedom. Under these circumstances, the CL of XdataThy ^ meaningful indicator of 

the level of agreement of the data and the theory. 

To take a simple example, suppose that there are 40 data and that Xdata wa — 30 
and XwAThy ~ ^o that, on the assumption of uncorrelated Gaussian errors, Eqn(2.6) 
gives Xdata,Thy = 46. The corresponding confidence levels are: xLa,M/A/d-o.f. = 30/39, 
CL= 0.849 ; x'wA,Thy/d.o.i. = 16/1, CL= 6.3 x 10"^ xLta,Thy/d.o.i. = 46/40, CL= 0.28. 
Thus the effect of the four standard deviation discrepancy observed in XwA,Thy i^ diluted 
to give an innocuous CL of 0.28 for the statistical estimator Xdata,Thy 

To take properly into account both the internal consistency of different measurements 
of the same quantity, and the level of agreement of the data with theory, a useful statistical 
procedure is to combine the confidence levels of the appropriate x^ functions. Since 
Xdata,WA X\vA,Thy independent statistical estimators, the corresponding CLs may 
be combined by use of the formula [21]: 

CL(q;i, — CKia2[l — \n{a-ia2)] (2.8) 

where ai and a2 are the two independent CLs to be combined. It follows that in the 
simple example considered above the combined CL has the value 5.8 x 10~^ so that the 
data/theory discrepancy is still well in evidence. Note that the combined CL is a factor 
493 smaller than the CL of XdataThy i^ ^^is case! In the following the combined CL given 
by Eqn(2.8) will be used to calculate the overall confidence level that the relevant data 
are consistent and that the data are in agreement with the SM, for different values of uih- 

In the above example each datum has the same sensitivity to the parameters of the 
theory. However, among the ^ 20 observables included in the global electroweak fits 
performed by the EWWG and the ~ 40 in the similar EWPDG fits, the majority are 
only weakly sensitive to the values of mn and This effect dilutes even further the 
hypothesis testing power of the the statistical estimator Xdata,Thy beyond that due to the 
dominant contribution of Xdata.wA discussed above. In the statistical analysis presented 
below, the separate contributions of Xdata wa XwA Thy ^o Xdata Thy "^iH tic extracted to 
provide separate answers to the questions posed in points (i) and (ii) above. The overall 
CL will then be calculated according to Eqn(2.8) above. In the case of a small number of 
sensitive observables ^it will be found that, unlike in the example discussed above, good 

^Indeed, for mn, there are only two such observables Ai and mw as discussed in Section 4 below. 
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agreement is found between the CL of the total x^: Xdata,WA + XwA,Thy combined 
CL calculated using Eqn(2.8). The former CL is then used as a statistical estimator for 
the indirect Higgs mass analysis. 

Previous authors [22, 23] have calculated normalised probability density functions 
(PDFs) giving the relative probability of different values of niH, by combining direct and 
indirect limits. Instead, in the present paper, the combined CL is found by combining 
the CLs of the direct and indirect measurements in region of overlap using Eqn(2.8). 
This combined CL gives an absolute rather than a relative probability that the SM is 
consistent with the data for any value of rriH- In this way the hypothesis testing aspect of 
the comparison of the data with the SM is addressed. This is not done by the normalised 
PDFs derived in References [22, 23]. 



3 Heavy Quark Asymmetry Measurements 



A discussion of the consistency of the b-quark asymmetry measurements in the data up 
to 1999 may be found in Reference [24]. The current LEP and SLD heavy flavour asymme- 
try measurements are collected in Table 1 (b-quarks) [25] and Table 2 (c-quarks) [25]. In 
Table 1 are reported eight independent measurements of the forward/backward asymme- 
try A^p^ as well as the direct SLD measurement of A^ from the forward/backward, left/right 
asymmetry. Table 2 contains seven LEP measurements of A^jf^ and the direct A,, mea- 
surement from SLD. For each LEP asymmetry measurement the corresponding value of 
Af, or Ac is estimated using the relation: 

= ^ iQ-b,c) (3.1) 

where Ai is the LEP+SLD average value of the charged lepton asymmetry parameter 
extracted by assuming charged lepton universality^: 



2viai 2ri 



vf + af 1 + ff 



(3.2) 



where 



The value used is [4]^ 



n = ^ = l-4sin'e|ff* (3.3) 
ai 



A = 0.1501(16) (3.4) 



The values of Ah and Ac derived in this manner are presented in the last columns of Tables 
1 and 2. The SM predictions for the values of A^ and Ac are 0.935 and 0.668 respectively. 



^The notation follows that of Reference [17] 

^Errors are quoted on the least significant digits, e.g. 4.123(32) means 4.123 ± 0.032. When two 
errors are quoted, the first is statistical and the second systematic 
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Experiment 


^FB 




ALEPH leptons 
DELPHI leptons 
L3 leptons 
OPAL leptons 
ALEPH inclusive 
DELPHI inclusive 
L3 jet-charge 
OPAL inclusive 
SLD 


0.1009(38)(17) 
0.1031(51)(24) 
0.1007(60)(35) 
0.0983(38)(18) 
0.1015(25)(12) 
0.0984(30)(15) 
0.0954(101)(56) 
0.1000(34)(18) 


0.896(39) 
0.916(51) 

0.895(70) 
0.873(38) 
0.902(27) 
0.874(32) 
0.847(103) 
0.888(35) 
0.925(14)(14) 


LEP average 
LEP-FSLD average 


0.0997(14)(7) 


0.885(14)(10) 
0.902(13) 



Table 1: The LEP and SLD measurements of b-quark asymmetry parameters. When two 
uncertainties are quoted, the first is statistical, the second systematic. 



Experiment 


^FB 


Ac 


ALEPH leptons 
DELPHI leptons 
L3 leptons 
OPAL leptons 
ALEPH D* 
DELPHI D* 
OPAL D* 
SLD 


0.0733(53) (36) 
0.0724(84) (62) 

0.0832(301)(197) 
0.0642(51)(37) 
0.0696(85) (33) 
0.0693(87) (27) 

0.0759(109) (57) 


0.651(57) 

0.643(93) 
0.739(320) 
0.570(56) 
0.618(81) 
0.615(81) 
0.674(109) 
0.670(20)(16) 


LEP average 
LEP-FSLD average 


0.0706(31)(17) 


0.627(26)(19) 
0.653(20) 



Table 2: The LEP and SLD measurements of c-quark asymmetry parameters. When two 
uncertainties are quoted, the first is statistical, the second systematic. 

with a negligible dependence on m,H and m,t at the scale of the present experimental 
errors. Also shown in Tables 1 and 2 are the LEP average values of A'p^, Af,, A'p^ and A^ 
as well as the LEP-I-SLD combined values of Af, and Ac. The uncertainties on the LEP 
average values of A^ and Ac come mainly from those on A^J^ and (1.6 % and 5.0 %) 
rather than that on Ai (1.1 %). The statistical and systematic errors on both the LEP 
average and the SLD measurements of Af, and Ac are of comparable magnitude. 

Values of the different estimators: Xdata,wA, XwA,Thy, and Xdata,Thy introduced 
above for the quantities Q = Ah, Ac are presented in Table 3. The XdataWA C'Ls of 0.92 
and 0.90 for A^ and Ac indicate good internal consistency of the data, but also, possibly, 
an over-estimate of systematic errors. The XwA,Thy C'Ls of 1.11 x 10"^ and 0.45 for A^ 
and Ac indicate in the former case a 2.5a discrepancy, and, in the latter, good agreement 
with the SM prediction. As in the example discussed above, the Ab discrepancy is not 
evident in the CL of XdataThy^ which takes the value 0.21. The combined CLs, according 
to Eqn(2.8), that the At, and Ac data arc both consistent and in agreement with the SM 
are 0.057 and 0.77 repectively. This would seem to indicate that it is not unreasonable 
that the deviation of Ai, from the SM prediction could be due to statistical fluctuation 
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xLa,wA/d-o.f., CL 




xLa,Tfev/d-o.f., CL 


Comb. CL 


A, 
Ac 


3.2/8, 0.92 
3.2/7, 0.90 


6.4/1, 0.011 
0.56/1, 0.45 


12.0/9, 0.21 
4.1/8, 0.85 


0.057 
0.77 



Table 3: Different estimators and CLs derived from the LEP and SLD measurements 
of At, and Ac. 



perhaps in combination with some unknown systematic effect. However, this conclusion 
requires confidence in the estimation of the systematic errors. As will be discussed below, 
there is some evidence, from the data itself, that the uncorrelated systematic errors may 
be somewhat overestimated, thus reducing from its true value the significance of the 
observed Aj, deviation. 

It is interesting to note that a goodness-of-fit estimator, indcpcndant of the test, 
is provided by the so-called 'Run Test' [26]. For the case of the At, data in Table 1, since 
all 9 independent measurements constitute a single 'run' (they are all less than the SM 
prediction) the corresponding CL is easily calculated. Since a single run can occur in only 
two ways (all data higher than or all data lower than the theoretical prediction) the CL 
is 2/2^ = 3.9 X 10~^. Unlike for the test, the CL of the Run Test is insensitive to over- 
or under-estimation of uncorrelated systematic errors. Since the CLs of the Run Test, 
of Xdata,WA of XwA,Thy independent, they may be combined into a single CL 

using the formula that generalises Eqn(2.8) to the case of three independent CLs: cti, a2 
and CKs [21]: 

CL(q;i, a2, as) = aia20(2[l - \n{aia2a2) + [^^("'i^^Qa)] j 

The combined CL for the Ai, data given by Eqn(3.5) is 2.5 x 10^^. The single run of the 
Ab data may be associated with a genuine deviation of the data from the SM prediction or 
a large correlated systematic effect of unknown origin. It is argued below that the latter 
explanation is unlikely. The third possible explanation, a statistical fluctuation, is also 
unlikely, given the small value of the combined confidence level. 

As there are ~ 10 independent measurements of both Af, and Ac, it is possible to 
compare errors estimated directly from the data, with the calculated statistical and esti- 
mated uncorrelated systematic errors on the weighted average values of Ai, and Ac shown 
in Tables 1 and 2 . The estimators for the error on the weighted average, a and its RMS 
uncertainty cr^ are given by the formulae [27]: 





-Q) 


\ N{N 


-1) 


a 





(3.6) 

<Ja = ; " (3.7) 

yJmN - 1) 

These formulae yield values of a of 0.0089(22) for A^, and 0.018(5) for A^ to be compared 
with the estimated total errors on the WA values of 0.013 and 0.020 respectively in these 
quantities. The agreement is good for Ac , but for A^ there is an indication, at the 1.9a 
level, that the uncorrelated systematic errors may be overestimated. This is confirmed by 



9 



calculation of the WA statistical error on the LEP+SLD weighted average value of Ab, 
which is just 0.0088, in perfect agreement with the value of a estimated directly from the 
data. Using this error to calculate XwA,Thy gives x^/d.o.f. = 14.0/1, CL = 1.8 x 10~^ a 
3.7(7 effect. Thus, although (see Tabic 1) the estimated systematic error on the LEP+SLD 
average value of is by no means the dominant one, the significance of the apparent 
deviation of At, from the SM prediction is very sensitive to it. There is some evidence, from 
the data itself, that the uncorrelated part of this systematic error may be over-estimated. 

The dominant source of correlated systematic error on both A^/^ and arises from 
the QCD corrections [28] of (2.96 ± 0.40)% for and (3.57 ± 0.76)% for A%. This 
source contributes 35% of the total systematic error on the LEP average value of A^, 
the remaining part being essentially uncorrelated between the different measurements. 
Thus the observed fractional discrepancy, 0.053, between the LEP average value of Ab 
and the SM prediction, is about 1.8 times larger than the QCD correction and about 13 
times larger than the estimated uncertainty on this correction. The latter would have to 
have been underestimated by more than an order of magnitude in order to explain the 
observed discrepancy with the SM prediction. This seems unlikely. Note, however, that 
the estimate of systematic error from the data itself using Eqns(3.6) and (3.7) gives no 
information on such a correlated uncertainty. 

In conclusion, the different measurements of A], are in very good agreement with 
each other, but their average value shows a -2.5(T deviation from the SM prediction. 
There is some evidence, from the data itself, that uncorrelated systematic errors may be 
over-estimated, thus possibly reducing the observed deviation from its true value. The 
correlated systematic error must have been underestimated by a large factor if the origin 
of the Ab deviation is unknown systematics rather than a breakdown of the SM. The 
measurements of Ac, on the other hand are found to be both consistent and in good 
agreement (within their much larger errors) with the SM prediction. Also however, as 
will be seen below, all of the hadronic asymmetries show similar fractional deviations from 
the SM parameters favoured by the purely leptonic data. The possibhty of deviations from 
the SM in the c-quark and light quark sectors as large as that observed in the b-quark 
sector is therefore not excluded by the asymmetry data. 

4 Sensitivities of Electroweak Observables to rrit and 
ruH 

To justify the restricted choice of observables used below to calculate the estimators 
for the data/SM comparison this Section presents some results of fits to obtain mn, or mn 
and nit. The overall approach used is the 'model-independent' one of References [16, 17, 
18] All charged lepton and heavy quark measurements from LEP and SLD are combined 
to obtain the independent observables: Ai, si, Ab, , Ac and Sc- The Af (/ = l,b,c) 
parameters are defined as in Eqn(3.2) above, with small additional correction terms in the 
case of Ab. The quantity, Sf, is defined as sj = v'j + a'j and so is proportional to the partial 
width for Z ^ ff decays. Again, due to the large mass of the b-quark, small corrections 
are included in this case [17]. The LEP-I-SLD average values of these observables used 
in the fits presented below are shown in Table 4. To take properly into account error 
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leptons 


c quarks 


b quarks 




Ai 


Sl 


Ac 




A, 


Sh 


Mcas. 


0.1501(16) 


0.25268(26) 


0.653(20) 


0.2897(50) 


0.902(13) 


0.3663(13) 


SM 


0.1467 


0.25272 


0.6677 


0.2882 


0.9347 


0.3647 


Dev(a) 


2.1 


-0.15 


-0.74 


0.3 


-2.5 


1.2 



Table 4: Measured values of Af and Sf compared to SM predictions for rut — 174 GeV, 
ruH = 100 GeV. Dev(a) = (Meas.-SM) /Error. 



X 




o-x 


{AX/ax)mt 


{AX/ax)m^ 


Ai 


0.1501 


0.0016 


3.1 


-5.1 


Sl 


0.25268 


0.00026 


2.4 


-0.11 


Ac 


0.653 


0.020 


0.10 


-0.18 


Sc 


0.2897 


0.0050 


0.18 


-0.077 


Ab 


0.902 


0.013 


0.012 


0.051 


Sh 


0.3663 


0.0013 


-0.12 


-0.22 


Su 


0.5014 


0.0015 


0.77 


-0.55 




1.3211 


0.0043 


0.93 


0.32 


mw 


80.426 


0.034 


3.5 


-3.1 



Table 5: Sensitivities of different measured quantites to rrit and tuh (see text). 

correlations, the directly measured values; Aj, — 0.925(20) and Ac — 0.670(26) from SLD 
are assigned speparate terms in the estimator. Correlations between Ai, and Ac 
resulting from Eqn(3.1) are included in the error matrix. Also shown in Tabic 4 are the 
SM predictions for rrit = 174 GeV and itlh = 100 GeV as well as normalised deviations. 
This Table has the same format and SM predictions as Table 3 of [17], with which it may 
be directly compared. 

The sensitivity of different observables to rrit and rriH is presented in Table 5 To 
take into account both the intrinsic sensitivity and the effect of experimental uncertainty, 
the quantities {AX/ax)mt ^-nd {AX/ax)m,H shown for each observable, X, with ex- 
perimental uncertainty ax- The quantity AX in {AX/ax)mt change in the value 
of X for a variation of rrit from 164 GeV to 184 GeV, with mn = 120 GeV and AX in 
[AX/(Tx)mH is the change in the value of X for a variation of mn from 20 GcV to 200 
GeV, with rrit = 174.3 GeV. Most of the sensitivity to rrit resides in the observables Ai, 
Sl and mw, to rrin in Ai and mw only. The sensitivity of si (or F;) to m^, but not to mn, 
has also been noted in a recent paper [29] . 

Also included in Table 5 are the observables: = vl + al and s'^j, (to be discussed 
below) which is similarly defined to Sc and Sb but for non-b quarks. Both these observables 
have a moderate sensitivity to both and mn- 

As pointed out above, if A^p^ is used as observable to estimate, via quantum correc- 
tions, mn, a very different value is obtained from that favoured by Ai or mw- This is due 
to the linear dependence of A^p^ on Ai (sec Eqn(3.1) above) and the 2.5a deviation of Af, 

''Note that the {AX/ax)mt entries of s; and Sb of the similar table in Reference [17] are incorrect 



11 



Source 


leptons 


b-quarks 


c- quarks 


r^had 
^FB 


hadronic mean 


overall mean 


Ai 


0.1501(16) 


0.1422(23) 


0.1423(72) 


0.1401(95) 


0.1421(21) 


0.1472(13) 



Table 6: Different experimental determinations, derived assuming tie correctness of tfie 
SM, of tfie leptonic asymmetry parameter Ai . 



Fitted (quantities 


iiiu [GeV] 


xVd.o.f. CL 


Ai{lept), si,Ab, Sb, Ac, Sc, mw 
Ai{all), m\.Y 

Ai{lept), si,Ab, Sb, Ac, Sc, mw mw'(NuTeV) 
Ai{aU), mw, mvi/(NuTeV) 
Ai{lept), mw 

Ai{lept), mw, mw^(NuTeV) 
Ai{had), mw 

Ai{had), mw, miy(NuTeV) 


y '-24 

97li 

113ti 

66^1^ 
1541^^ 

19611 


14.7/8, 0.065 

1.99/1, 0.16 
23.3/9, 0.0056 
10.6./2, 0.0050 

0.20/1, 0.66 
10.9/2, 0.0043 

8.0/1, 0.0047 
14.6/2, 0.00068 



Table 7: Results of mn fits to different sets of observables. 



from the SM prediction discussed in the previous Section^. Assuming the correctness of 
the SM, 'hadronic' values of Ai may be extracted from the measurements of A^p^ and A^jf^ 
by substituting the SM predictions for A^ ,Ac (which are essentially independent of mn 
and mt) into Eqn(3.1). Another, independent 'hadronic' value of Ai may be derived from 
the value of sin^ ^t// obtained from the SM analysis the quark anti-quark charge asym- 
metry, [4]- These different 'hadronic' determinations of Ai, obtained by assuming 
the correctness of the SM, are presented, together with the 'leptonic' value from Table 4, 
in Table 6. Note that the 'leptonic' value of Ai, although derived assuming charged lepton 
universality does not assume the correctness of the SM, only that the process Z ^ U is 
described by some effective vector and axial vector couphngs so that Eqn(3.1) is obeyed. 

The 'leptonic' value of Ai in Table 6 is derived from e, n and r forward/backward 
charge asymmetries and from r-polarisation measurements. The hadronic ones from 
quark forward/backward charge asymmetries. In fact the Ai derived uniquely from r- 
polarisation measurements: Ai{t — pain) = 0.1465(33) lies almost exactly mid- way be- 
tween the SLD ALR and LEP weighted average of 0.1513(19) and the value of 
Aiijiad) quoted in Table 6. It is 1.3a below the former and l.lo" above the latter. In the 
1996 data set [30] the difference between Ai{t — pain) and the ALR, J^p^ average was 

*This effect is particularly transparent in Fig.l of Reference [18] or in Fig 15.1 of Reference [4] 





all data except NuTeV 


all data 


This 
paper 


rriH 
mt 

xVd-O-f, CL 


102li 

175.0it2 
14.7/8, 0.065 


107±1« 
173.7ltl 
23.3/9, 0.0056 


EWWG 
Ref. [4] 


rriH 

mt 

xVd-o.f, CL 


^-■--36 

175.3;t^ 
16.7/14, 0.27 


174.3lti 
25.4/15, 0.045 



Table 8: Global electroweak fits for mn and mt- 
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much larger, 2.5a, so that the inclusion (or not) of the r-polarisation data had a large 
effect on the value of Ai extracted using Eqn(3.1). This was discussed in some detail in 
Reference [16]. In another paper discussing the same 1996 data set [31] it was pointed 
out that, considering also the r-poIarisation data as 'hadronic' (because of the predom- 
inantly hadronic final states), the value of sin^ Ogj'j derived from the leptonic ALR and 

A^p^ measurements was found to differ by more than 3a from that given by the 'hadronic' 
ones, i.e. r-polarisation and quark asymmetry measurements. The situation is much 
improved in the current (essentially final) LEP-I-SLC data. Excluding the r-polarisation 
measurements gives a minor change in the WA Aj, value: Ab{T — poln out) = 0.898(13) to 
be compared with the value quoted in the last row of Table 1. The deviation of Ab from 
the SM prediction is only increased from 2.5a to 2.8a, instead of the ~ la increase found 
in the 1996 data set [16]. 

Because of its strong dependence on ttih and nit then, unlike in the case of A^ and 
Ac, no definite SM prediction exists for Ai. However it is of interest to compare the 
'leptonic' value of Ai, Ai(lept) with the different 'hadronic' values Ai{had). The follow- 
ing values and confidence levels are obtained: XhadWA/^-^-^- — 0.047/2, CL= 0.977; 
xlu,wA/d.o.l. = 9.0/3, CL= 0.029: xLd,/ept/d.o.f. = 9.2/1, CL= 0.0024. Thus the three 
hadronic determinations are very consistent with each other, whereas the hadronic and 
leptonic determinations differ by 3 standard deviations. This poor overall consistency of 
the different values of Ai, extracted assuming the correctness of the SM for the quark 
couplings, must be taken into account when assessing the overall level of agreement of 
the data with the SM.^ It is important to stress that this mismatch is not the result of 
any inconsistency evident in the experimental data themselves, but rather the result of 
interpreting the data according to the SM prediction. 

The following strategy is now followed for fits to obtain hmits on mn- In a first 
step, fits similar to those previously presented in References [17, 18] arc performed to the 
the entire LEP+SLD data set contributing to the six observables of Table 4, as well as 
LEP-I-FERMILAB combined direct measurement of mv^: mw = 80.426(34) GeV. Other 
fits are done including also the indirect determination of mw' miy(NuTeV) = 80.136(83) 
GeV by NuTeV [32]. Only niH is varied in the fits, the other important parameters: mz, 
rrit, a{mz) and as{mz) being fixed at their measured values of 91.1875 GeV, 174.3 GeV, 

0. 007755 and 0.118 respectively. The effect of variation of the second and third of these 
parameters, within their experimental uncertainies, on the CL for agreement of the data 
with the SM, will be discussed in Section 5 below. The values of other fixed parameters 
are specified in References [16, 17, 18]. 

In the fits, a numerical parameterisation, accurate at the per mil level, of the effective 
weak mixing angle given the two- loop ZFITTER 5.10 program [33] was used^°: 

sin^ Ogjj = 0.233657 - 8.42 X lO'^m^ - 3.86 x 10"^ In rrit + 5.00 x IQ-^lnrriH (4.1) 

where rrit and rriH are expressed in GeV units. The overall normalisation factors Pf {f = 

1, V, u, d, h) for fermionic widths of the Z are given by a numerical parametrisation similar 

^Indeed, the consistency of the three different 'hadronic' estimates of Ai is much better than expected. 
Because of the large statistical uncertainties of the b- and c-quark data this is most likely due to a chance 
co-incidence rather than any over-estimation of systematic errors. 

^°The formula is valid at the quoted accuracy for wih > 40 GeV. For lower values of niH, small 
corrections are made to the constant term and the coefficient of In niH ■ 



13 



to Eqn(4.1) of the entries in Table 2 of Reference [34]. For mw, the parameterisation of 
Reference [35] was used. 

The fits for niH arc then repeated using only the 'mi^-scnsitive' observables Ai{all) 
and niw where Ai{all) is the weighted average of the leptonic and hadronic values of Ai 
given in the last column of Table 6. Similar fits are performed including also mvF(NuTeV). 
The results of this comparison are shown in the first four rows of Table 7. In can be seen 
that essentially the same range of Higgs masses is obtained whether fits are made to the 
complete set of electroweak observables or only to the m/j-sensitive ones Ai and rriw- The 
fit results presented in the fifth and sixth rows of Table 7 demonstrate that very low values 
of niH, with best fit values incompatible with the 95% direct lower hmit of 114.4 GeV, are 
found when fitting only the m//-sensitive observables A{lept)i and niw As shown in the 
last two rows of Table 7, much higher values of mn are found when A{lept)i is replaced 
by A{had)i in the fits. This is a consequence of the deviation of the measured value of Ab 
from the SM expectation, and the strong correlation between Ai and Ab resulting from 
Eqn(3.1), when A^p^ is measured. In all cases inclusion of the NuTeV mw measurement 
results in slightly higher fitted values for niH and reduces all confidence levels, by about 
an order of magnitude, to values less than 0.01. 

As another cross-check of both the fitting procedure and the uih, rrit sensitivity of 
different observables, simultaneous fits to ruH and mt were performed including also in 
the estimator the directly measured value of rrit from FERMILAB: mt = 174.3(5.1) 
GeV. The results of these fits, both including and excluding the indirect NuTeV m^ 
measurement, are presented in Table 8, where they may be compared with the results 
of similar fits from the most recent EWWG report [4]. Slightly lower fitted values of 
mn are found in the latter, probably due to the inclusion of other observables such as 
Fz, and Qw{Cs) from atomic parity- violating experiments, that have some sensitivity to 
mfj, in the EWWG fits. The uncertainties on both mn and mt found in the two sets 
of fits are very similar. In fact slightly more precise values of mH are obtained in the 
fits of the present paper. This, in combination with the results shown in the first four 
rows of Table 7, shows that the restriction to Ai and mw entails no significant loss of 
sensitivity in the indirect determination of mn- The dilution effect discussed in Section 
1 of the hypothesis-testing power of the estimator, due to the inclusion of unaveraged 
equivalent obervables, or additional 'noise' observables, that are insensitive to mn and 
mt, is evident in the x^/d.o.fs and CLs of the fits that are also presented in Table 8. 
The EWWG fits have a CL that is a factor of 8(4) times larger than those of the present 
paper for the fits including (excluding) the NuTeV mw- measurement. A more detailed 
discussion of the global EWWG fits is found in Section 8 below. 

5 Combining Confidence Levels of Direct and Indi- 
rect Limits on itlh 

The combined result of the direct searches for the Standard Model Higgs Boson by the 
LEP Collaborations ALEPH, DELPHI, L3 and OPAL is given in Fig.9 of Reference [10]. 
This shows the confidence level ratio: CLg = CLg+b/ CLb as a function of mn- C^s+b is the 
confidence level of the signal-plus-background hypothesis and CL^ that of the background- 
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Fi gure 1: Illustration of the conibina.tioii of direct (Cljs+b) sud indirect CL()(^g]y[ + Xwa) 
tuh confidence levels using Eqn(2.8). The observables used to calculate CL(xsm + Xwa) 
are Ai{aU) and mw- 
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Figure 2: Combined nin confidence levels. The observables used to calculate 
CL(XsM + Xwa) are Ai and mw 
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Figure 3: Combined uih confidence levels. The observables used to calculate 
CL(xiM + Xwa) are A, mw Bind mvK(NuTeV). 



Observables 


mn [GeV] 


120 


160 


200 


240 


280 


Ai{all), mw 
x'wA/d-o.f. = 9.0/3 
CL = 0.029 




0.51 


0.35 


0.13 


0.047 


0.0074 


CHx'sM + Xwa) 


0.065 


0.049 


0.022 


0.0073 


0.0021 


CL{Comb) 


0.076 


0.057 


0.024 


0.0104 


0.0020 


Ai{lept), mw 
xlvA/d.oI. = 1.6/2 
CL = 0.45 


CL(xIm) 


0.30 


0.054 


0.008 


0.0011 


0.00017 


CUxsM + Xwa) 


0.41 


0.11 


0.024 


0.0044 


0.00079 


CL{Comb) 


0.41 


0.11 


0.024 


0.0044 


0.00079 


Ai{had), mw 
X^^/d.o.f. = 0.047/2 
CL = 0.98 


cHxIm) 


0.11 


0.024 


0.029 


0.026 


0.019 


CHXsM + Xwa) 


0.058 


0.111 


0.130 


0.118 


0.091 


CL{Comb) 


0.058 


0.111 


0.130 


0.118 


0.091 



Tabic 9: Confidence levels as a function of mn for different sets of observables. The 
values of x^^/d.o.f. and CL refer to Ai. 
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Observables 


rriH [GeV] 

-'^ L J 


120 


160 


200 


240 


280 


Ai(all). Tfiu/ . 
mvi/(NuTeV) 


VA5M/ 


0.012 


0.015 


0.0089 


0.0036 


0.0012 


CL(y^,^ + Ywa) 


0.0027 


0.0033 


0.0021 


0.00096 


0.00036 


CL(ro/;/M 


0.0030 


0.0037 


0.0021 


0.00107 


0.00039 


mvi/(NuTeV) 




0.0071 


0025 


00066 


00013 

\J iXJ \J\J -L ^ 


3.3 xlO"^ 


CL (X5M + Xwa) 


0.018 


0.0072 


0.0022 


0.00056 


0.00014 


CL{Comh) 


0.022 


0.0089 


0.0027 


0.00065 


0.00018 


A/(/ia(i), mw, 
mw (NnTeV) 




0.00031 


0.0012 


0.0022 


0.0029 


0.0029 


CL(x'9M + xivA) 


0.0021 


0.007 


0.012 


0.015 


0.015 


CL(C'o/;/6) 


0.0028 


0.0091 


0.016 


0.020 


0.019 



Table 10: Confidence levels as a function of uih for different sets of observables. 





-4/(.4//) 


Aiilcpt) 


Aiiliad) 


niH (GcV) 








111 


8.4 X 10"^ 


6.3 X 10-6 


6.6 X 10-^ 


113 


3.2 X 10-^ 


2.1 X 10-3 


2.6 X 10-^ 


115 


0.047 


0.23 


0.041 


140 


0.062 


0.23 


0.088 


180 


0.034 


0.054 


0.13 


220 


0.013 


0.010 


0.13 


260 


3.9 X 10-3 


1.9 X 10-3 


0.11 


300 


1.1 X 10-3 


3.3 X 10-^ 


0.077 



Tabic 11: Combined confidence levels CL for consistency witli tlie SM as a function of 
tuh- Observables used in the estimator: Ai and mw- 





Ai{All) 


Ai{lept) 


Ai{had) 


mn (GeV) 








111 


3.6 X 10-s 


3.0 X 10-^ 


2.6 X 10-8 


113 


1.6 X 10-5 


1.1 X 10-^ 


1.1 X 10-5 


115 


2.8 X 10-3 


0.017 


2.2 X 10-3 


140 


3.3 X 10-3 


0.012 


4.2 X 10-3 


180 


2.8 X 10-3 


4.0 X 10-3 


9.8 X 10-3 


220 


1.5 X 10-3 


1.1 X 10-3 


0.014 


260 


6.0 X 10-^ 


2.8 X 10-^ 


0.015 


300 


2.1 X 10-^ 


6.6 X 10-5 


0.014 



Table 12: Combined confidence levels CL for consistency with the SM as a function of 
mn- Observables used in the estimator: Ai, mw cLnd m^(NuTeV). 
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only hypothesis. Inspection of the figure shows that, at percent level accuracy: CL^ = 
10-^,0.05,0.08 for mn = 111,114.4,120 GeV, respectively. For the present study it is 
preferred to work directly with CL^.^^, which is similar to the confidence level given 
by comparing the SM to Z-decay data, to obtain indirect rriH limits. As shown in Fig. 7 
of Reference [10], the value of CL^ is about 0.8 in the region 110 GcV < uih < 120 GeV, 
of interest for the present study. This gives the estimates: CLg+b = 10~^, 0.04, 0.64 for 
ruH = 111, 114.4, 120GeV, respectively. 

The following numerical parameterisation of CLs+b is used: 

111 GeV < ruH < 114.4 GeV 
logCL,+5 = 1.382mjy(GeV) - 159.51 (5.1) 
114.4 GeV < ruH < 120 GeV 

As shown in Fig.l, the function of Eqn(5.2) has the same value and first derivative as 
that of Eqn(5.1), at the matching point itlh = llAA GeV, and vanishing first derivative 
at niH = 120 GeV, where CLg+b = 0.64. Allowing for the overall scale factor of 0.8, 
the parameterisation of Eqns(5.1) and (5.2) describes well the experimentally determined 
curve of CLs+ft/CLf, in Fig.9 of Reference [10]. It should be noted, however, that the 
precise shape of CLs+;, has only a small effect on the final confidence level curves to be 
presented below. The direct search excludes, with a CL of < 10~^, the possibility that 
the SM Higgs boson exists with mass of less than 113 GeV, and gives essentially no 
information for mn > 115 GeV. Thus the region where it is of interest to combine CLs+b 
with indirect confidence levels covers only a narrow range of mn- 

In order to define the confidence level for agreement of Z-decay data with the SM, 
the of the data/SM comparision is simply calculated as a function of uih, setting rrit 
and a{'mz) to the measured values given above. Therefore no fit to the data is necessary. 
The sensitivity of the CL curves to the assumed values of rrit and a{mz) is discussed 
below. In this Section only the 'mi^-sensitivc' obscrvablcs Ai and mw are included in 
the estimator, where the W mass is either the directly measured value from LEP and 
FERMILAB or the indirectly determined NuTeV value. Ai is determined either by using 
all asymmety data {Ai(all)), lepton data only {Ai{lept)) or only hadronic data {Ai{had)). 
The corresponding values arc presented in Table 6 above. To take into account the internal 
consistency of the different data sets the values of Xaii,WAy xlept,WAi x\ad,WA ^-re added 
to the x^ of the SM comparison: Xw^a sm Q^ich. case. As shown below, almost identical 
CLs are found using either Xx,wa + Xwa,sm — lept, had) or by combining the 
CLs of Xx,wa Xwa,sm using Eqn(2.8). The former CL is then combined with CLs+6 
using Eqn(2.8) to yield the direct plus indirect confidence level curves shown below. The 
values of x^/d.o.f. for XhadWA ^-nd xliiwA given above; that for the Icptonic data: 
xfept,wA/^-^-^- — 1-6/2, CL = 0.45, given by combining the Ai values obtained from 
lepton forward/backward asymmetries and tau polarisation measurements from LEP and 
the Alr measurement from SLD, is taken from Reference [4]. 

Some typical CLs for the indirect rriH analysis, obtained as described above, are pre- 
sented in Table 9 (observables considered: Ai, mw) and Table 10 (observables considered: 
Ai, mw, mvF(NuTeV)). In all cases good agreement is found between CL{x%m~^Xw a) 
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CL{Comb) calculated using Eqn(2.8), where the abbreviations Xsm = Xwa,sm Xwa = 
Xjc wa have been introduced. 

The combination of G\j{x\m ^ Xw a) (indirect measurements) and CLs+b (direct mea- 
surements) for different values of tuh is illustrated in Fig.l. Since CL^+j, provides lit- 
tle information on vriH for mn > 114.4 GeV (the 95% CL lower limit of the direct 
search), GIj{x%m + Xwa) combined with Ghg+b provided that CLiComh) is less than 
CIj{x%m ^Xwa)- the contrary case G\j{x%m ^Xwa) alone is used. Thus the algorithm 
used to obtain the combined confidence level, CL, is: 

• Calculate 0:3 from cci = CL{x^y[ -\- x^a) 0:2 = CL5+6 according to Eqn(2.8). 

• If as < CL(x|a^ + Xwa)^ set CL = CL(ai aa). 

• If as > CL(x|m + Xwa), set CL = CL(x|m + Xwa)- 

In Fig.l the dotted curve shows CLs+6, the dashed curve CL(x|jy^ -|- Xwa) the solid 
curve CL. 

Curves of CL calculated in this manner are shown in Figs. 2 and 3. When Ai{lept) 
and ^4^(0//) are used, the general shape, with a sharp peak above, but close to, the direct 
lower limit and a rapid fall-off for higher values of itih is similar to that of the PDFs 
presented in References [22, 23]. However CL, unlike the PDFs, gives an estimate of the 
absolute probability that the data is consistent with the SM for a given value of mn- 
An exception to this behaviour is provided by the data sets Ai{had), mw and Ai{had), 
mw, 7TT'W'(NuTeV) where the maximum of CL occurs at much higher values of uih and 
on average much larger values of CL are obtained. The confidence level curves in Figs. 2 
and 3 are presented in numerical form in Tables 11 and 12, in the same format as the 
similar curves considered in Section 7 below. The latter use all precision observables , 
rather than only the mif-sensitive ones as in Figs. 2 and 3. Comparison of the two sets of 
curves then shows the effect of 'non mi^-sensitive' observables on the level of agreement 
with the SM prediction. 

The effect on the CL curves of variation of the value of rrit by plus or minus the 
experimental error around the measured value is shown in Fig.4 and Table 13. The effect 
of a similar variation of a{mz) is shown in Fig. 5 and Table 13. For large values of rriH 
this variation of changes the values of CL by many orders of magnitude. Even so, in 
the case of the data set Ai{lept), mw, shown in Figs. 4 and 5, which can be argued (see 
below) to be likely to give the most rehable estimate oi uih: CL is still only, at best, 
0.01, at rriH = 300 GeV. Change of a{mz) within the current experimental errors can 
change CL by up to an order of magnitude, but the effect is much less dramatic than for 
nit. Clearly a much improved measurement of is needed to significantly improve the 
indirect limits on tjih- 

The analysis presented in this Section has many similarities with that of Reference [19] 
where confidence levels taking into account both direct and indirect information on itlh 
were derived. However, the present writer has doubts about the mathematical correctness 
of the method used in Reference [19]. It consists of combining the confidence level given 
by the Xmin d.o.f. of a fit with that derived from Ax^ = x'^i'n^H — 114 GeV) —Xmin 
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Figure 4: Dependence of combined mu conRdence levels on the value of m^. The 
observables used to calculate CL(xgi^ + Xwa) Aiilept) and mw- ot{mz) = 0.007755 is 
assumed. 
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140 


0.032 


0.61 


0.46 


0.073 
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3.4 X 10-3 


0.30 


0.18 


0.010 


220 


3.5 X 10-^ 


0.11 


0.051 


1.3 X 10-3 
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3.6 X 10-5 


0.036 


0.013 


1.7 X 10-^ 
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4.0 X 10-6 


0.010 


3.1 X 10-=^ 


2.3 X 10-5 



Table 13: Combined confidence level curves CL for variation of rrit and a{mz) by plus 
or minus one standard deviation around their measured values. Observables used in the 
estimator: Ai{lept) and m^w 



20 



0.007752 

0.007755 

0.007758 



no 130 150 170 190 210 230 250 270 290 310 
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Figure 5: Dependence of combined mu confidence levels on the vsdue of a{mz)- The 
observables used to calculate CL(xg]y[ + x^^) are Ai{lept) and mw- iTit — 178 GeV is 
assumed. 

the same fit, assuming them to be independent. Certainly, the formula used to combine 
the confidence levels, simple multiplication instead of Eqn(2.8) above, is incorrect^^. The 
approach used is essentially to replace the measured CLg+b curve by a ^-function at 
m,H = 114 GeV. As can be seen in Fig.l, this is quite a good approximation. The global 
fits used in Reference [19] took no account of the dilution of the hypothesis testing power 
of the estimator resulting from the use of unaveraged and insensitive observables, as 
discussed in Section 2 above and, in more detail, in Section 8 below. For instance there 
is the statement: 'The global SM fit was excellent in 1998 and has now (2002) become 
poor'. This is not at all true of the contribution to the of the m//-sensitivc observables, 
which is similarly high for both data sets. In fact, the high confidence level for the 1998 
global fit is a consequence of an anomalously low contribution from the non m^f-sensitive 
observables [18]. Also, although the NuTeV measurement was discussed, together with 
Ab, as a possible source of anomaly relative to the SM prediction, only the published 
interpretation as a measure of sin^ Qon-sheii considered. Also the correctness of 

the SM prediction for the quark couplings to the Z was assumed to define different values 
of sin^ 0g'J_^ (equivalent to Ai in the present paper). The same two assumptions have 
been made in the analysis presented in this Section. In the analysis of model-independent 
observables presented in Section 7 below the alternative interpretation of the NuTeV 
result, discussed in the following Section, as a measurement of the Zuu coupling is also 
used, where no a priori assumptions are made concerning the couplings of the Z to fermion 
pairs. 

^^The present author made the same mistake in Reference [16] 
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An important point stressed in Reference [19], first pointed out in Reference [18], is 
that, regardless of how the Ai, anomaly is interpreted (statistical fluctuation, unknown sys- 
tematic effect or new physics) the most reliable estimate of niH must be that derived from 
the charged lepton asymmetries and rriyy that give consistent predictions for this quan- 
tity. Inclusion of the ^4^^ measurements in the fit results in a positive ~ 50 GcV bias 
on the 95% CL upper limit on niH, due to the Ai-Ab correlation resulting from Eqn(3.1). 
Although fits with and without the 'anomalous' NuTeV measurement are routinely pre- 
sented by the EWWG, fits excluding the (equally 'anomalous') hadron asymmetry data 
that, would provide the most reliable estimate of mj^, have (to my best knowledge) never 
been shown in the periodical updates of the status of clectroweak measurements, such 
as Reference [4], produced by this working group. I agree with almost all of the general 
conclusions of Reference [19], in particular, that if the lepton asymmetry data is correct 
then, regardless of the status of the hadronic asymmetry data, the SM provides only a 
poor global description of the data. 1 would also remark that the confidence levels of 
the global fits quoted, although small, will become even smaller when corrected for the 
dilution effects discussed in Section 2 above. Use of the correct formula for combining 
confidence levels will, on the other hand, give higher combined confidence levels. The 
analysis of Reference [19] gives, however, no hint of the very low values of CL for large 
values of TTiH apparent in Figs. 2 and 3 and Tables 11 and 12. Finally, in connection with 
Reference [19], as discussed in Section 8 below, the use of Ax^ to provide confidence 
levels for parameter estimation is of doubtful validity when, as is the case for the current 
electroweak data, the absolute confidence level derived from Xmin shows that the model 
containing the parameter of interest does not adequately describe the data. 



6 Alternative Interpretations of the NuTeV Experi- 
ment 

The pubhcation of the results of the NuTeV experiment [32] gives an estimation of 
the value of the on-shell weak mixing angle: 

sin^ e^-'''^^^ = 0.2277(13) (9) (6.1) 

that may be translated directly into a W-mass measurement via the defining relation of 
the on-shell renormalisation scheme: 

sin2^^^-'^^" = l-l4^ (6.2) 

This interpretation however requires that the Zi/V coupling is assigned its standard model 
value. As discussed in Reference [36] and shown, for example, in Fig.l of Reference [38], 
the experiment actually measures a quantity that is sensitive both to sin^ Qon-sheii 
the strength of the Zuu coupling, which may be specified by a parameter po such that: 

s^^vl + al^plMSM) (6.3) 
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Assuming the SM value (po = 1) gives a prediction for sin^^^ 

sin^ e^-''^"" = 0.22733(135)(93)-8.8xl0-^(mt[GeV]2-1752)+3.2xl0-^ln(mH[GeV]/150) 

(6.4) 

Assuming instead the SM value of sin^ 6^ the experiment provides a measurement 
of po [36]: 

po = 0.9942(13) (16) + 2.4 x lQ-\mt[GeV]'^ - 175^) - 1.6 x 10-Mn(m^^[GeV]/150) (6.5) 

Choosing the values rrit — 175 GeV and itih = 150 GeV consistent with the measured 
value of niw, (6.5) gives a measurement of the model-independent parameter Su'- 

s^(NuTeV) = 0.4992(21) (6.6) 

which may be compared with the LEP measurement quoted in Table 5: 

s^(LEP) = 0.5014(15) (6.7) 

Since: 

s^(LEP) - s^(NuTeV) = 0.0022(26) (6.8) 

the two measurements are quite consistent {XwAs^/^-'^-f- — 0.78/1, CL = 0.38) and yield 
the weighted average value: 

s^(LEP + NuTeV) = 0.5006(12) (6.9) 

Alternatively assuming that po = I and using Eqns(6.2) and (6.4) to obtain mw gives: 

s^(NuTeV) = s,.(SM) = 0.5050 (6.10) 
mH.(NuTeV) = 80.136(83) GeV (6.11) 

In this case the assumed, SM, value of Si, differs from the LEP measurement by 2Aa 
{x^/d.o.f. = 5.76/1, CL = 0.016) and also 

mvK(LEP + FERMILAB) - mv^(NuTeV) = 0.290(90) (6.12) 

the 3.2(7 discrepancy {x^/d.o.f. — 10.4/1, CL = 0.0013) mentioned in Reference [32]. 
Using Eqn(2.8) to combine the data consistency CLs for Su and mw yields an overall 
CL of 2.5 X 10^^. Thus on the assumption that the NuTeV measurement is correct, the 
alternative interpretation of the experiment is strongly favoured statistically as the ratio 
of data consistency CLs of the two interpretations is ~ 1.5 x 10^. 

For both interpretations the SM prediction is unfavoured. For the standard one {mw 
measurement and po = 1) the CL of the SM comparison is that just quoted: 2.5 x 10~^. 
For the alternative interpretation (po ^ 1) it is found that: 

s^(LEP + NuTeV) - s^(SM) = -0.0044(12) (6.13) 

or a 3.7(7 {x^/d.o.f. = 13.4/1, CL = 2.5 x 10"^) deviation from the SM prediction. The 
alternative interpretation thus shows exactly the same deviation from the SM as the one 
proposed in Reference [32]. 
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A number of theoretical assumptions must be made in order to derive Eqns(6.4) and 
(6.5) from the experimental quantities: 



a{uFe vX) 



VX) 



a{vFe^li-Xy 



a{vFe 11+ X) 



actually measured by the NuTeV experiment. A recent concise review of the situation 
may be found in Reference [20] in which citations of related work can be found. The most 
important and extensively discussed assumption concerns the supposed symmetry of the 
strange sea momentum distribution in a nucleon. A recent analysis by the CTEQ Col- 
laboration [37] presented in Reference [20] finds some evidence for a positive asymmetry 
of the strange quark sea: 



It is pointed out in Reference [20] that an asymmetry of 0.002 has the effect of reducing by 
42% the discrepancy between the measured value of sin^ Qon-sheii jgj^jyg^j from the direct 
mw measurements, and the value of the same quantity found from Eqn(6.4). The alterna- 
tive, and statistically favoured, interpretation of the NuTeV experiment as a measurement 
of po, was not considered in Reference [20], but in view of the linear correlation between 
sin^ and po provided by the measurement^^, it is reasonable to suppose that, in 

the alternative interpretation, the strange quark sea asymmetry will reduce the deviation 
of Po (or, equivalently s„) from the SM expectation by the same fraction. Thus the esti- 
mated value of Siy, correcting for the efi:ect of the asymmetry of Eqn(6.14) is 0.5015(21), 
which agrees prefectly with LEP measurement in Eqn(6.7). The LEP-|-NuTeV weighted 
average becomes 0.5014(12), which still hes 3.0(7 below the SM expectation. The appar- 
ent anomaly in the ZuV coupling is therefore reduced, but not removed, by the estimated 
effect of a strange quark sea asymmetry on the NuTeV results. 

It remains true however that because of the many systematic effects, detailed in Ref- 
erence [20], the results of the NuTeV experiments are less 'sure' than the measurement 
of the related quantity Fj^j, at LEP. Because of this some authors [43] prefer to adopt a 
conservative position and exclude the NuTeV results completely from global electroweak 
analyses. In contrast, the present paper takes a strictly neutral position on the question 
of the reliability, or otherwise, of the NuTeV results. In the following Section then, CLs 
as a function of niH will be calculated using all available LEP and SLD data on the as- 
sumption of either of the two possible interpretations of the NuTeV experiment, or by 
excluding the experiment. The CLs obtained are later compared in Section 8 below with 
those obtained from the global EWWG and EWPDG fits. 

7 Model-Independent Observables Compared to SM 
Predictions 

Following the model-independent approach of References [16, 17, 18] essentially all 
precision information on the SM provided, to date, by the LEP, SLC and FERMILAB 

^^See, for example, Figure 1 of Reference [38]. 




(6.14) 
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experimental programs, as well as by the NuTeV experiment, is contained in the values 
of the nine observables listed in Table 14. The other important quantities rrit and Q;(mz) 
are not included as they are here considered as input parameters for the SM prediction 
rather than measurements which provide a test of the SM. All information from leptonic 
forward/backward charge asymmetry and r-polarisation measurements as well as the SLC 
ALR measurement is condensed into the single parameter Ai{lept), equivalent to sin^ ^i//- 
The quantity si is derived from the width of the Z for decay into charged leptons on the 
assumption of charged-lepton universality. The quantities Ac and Ai, are derived using 
Eqn(3.1) from the c- and b-quark forward/backward charge asymmetries measured at 
LEP as well as from the SLC measurement of forward/backward-left /right asymmetries 
of c and b quarks. Sc and Sb are obtained from the Z decay widths into c and b quarks, 
more conventionally expressed in terms of the ratios: Rq — rQ/Fhad, {Q — c,b), using 
the relation: 

The QED and QCD correction factors are given in Reference [16]. The observable s,^ 
is given by the invisible width, Fj^^, of the Z boson, determined from the Z-boson total 
width, r^, the hadronic width, Thad and the leptonic width, F/ via the relation: 

^inv — ^Z — F/iad — 3F/ (7.2) 

Since the hadronic width of the Z is quite precisely measured: Thad = 1.7444(20) GeV, 
the measurements of Rb or Sb can be used, in combination with the former, to extract the 
quantity: 

q=u,c q=d,s 

The subscript 'nb' here stands for 'non-b' quarks. As will be discussed below, the mea- 
surements of Thad and s^b provide much more stringent constraints on the possible values 
of the couplings of non-b down- type quarks to the Z than the existing direct measurements 
of these couplings, which have large experimental errors. 

The experimental errors on the observables listed in Table 14 are largely uncorrelated 
between the observables, which facihtates calculation of a estimator for global SM 
comparisons. Correlations exist between: Ai, AcilMF) and ^^(LEP) due to the use of 
Eqn(3.1) to obtain ^^(LEP) and ^^(LEP). Because of the small uncertainty on T^ad the 
errors on Sb and s'^j, are strongly anticorrelated. Weaker correlations exist between Sc and 
Sb- In view of the relatively poor precision of the Sc measurement in comparison with 
those of Sb and s^^^, and the correlations between these three observables, the contribution 
of the former is omitted from the estimator used in the global comparisons with the 
SM shown below. As previously mentioned, to take properly into account correlations, 
the direct SLC measurements of A^. and A^ are assigned separate terms from the LEP 
measurements in the estimator. In Table 14, however, the weighted average LEP-I-SLC 
values of A^. and Ai, are quoted. 

The value of 'Sj^ given in Table 14 is the LEP+NuTeV weighted average, i.e. the 
'alternative' interpretation of the NuTeV experiment is taken. This is mandatory in a 
model-independent analysis where no a priori assumptions concerning the couplings of 
the leptons or quarks to the Z (with the exception of charged lepton universality, well 
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^expt/o'X (%) 






Alilept) 


0.1501(16) 


1.07 


0.1481 


1.05 


Sl 


0.25268(26) 


0.103 


0.25277 


-0.35 


Ac 


0.653(20) 


3.06 


0.668 


-1.06 




0.2897(50) 


1.73 


0.2884 


0.26 


At 


0.902(13) 


1.44 


0.9347 


-2.51 


Sh 


0.3663(13) 


0.35 


0.3648 


1.15 


Su 


0.5006(12) 


0.24 


0.5050 


-3.7 




1.3211(43) 


0.33 


1.3218 


-0.16 


mw 


80.426(34) 


0.042 


80.394 


0.94 



Table 14: Model-independent electroweak observables. The SM predictions correspond 
to mt = 178 GeV, mn = 120 GeV and a{mz) = 0.007755. 



Coupling 


Expt value 


SM 


(Exp-SM)/(T 


Vl 


-0.03783(41) 


-0.03734 


-1.2 


ai 


-0.50125(26) 


-0.50137 


0.46 




0.5003(6) 


0.50251 


-3.7 




0.1875(69) 


0.1921 


-0.67 


flc 


0.5045(50) 


0.5015 


0.60 


Vh 


-0.3232(78) 


-0.3435 


2.6 


Ob 


-0.5133(50) 


-0.4983 


-3.0 



Table 15: Effective vector and axial-vector coupling constants of the Z boson to lepton 
neutrino and heavy quark pairs. The SM predictions correspond to rrit = 1 78 GeV, uih 
= 120 GeV and a{mz) = 0.007755. 

respected by experiment) is made. The NuTeV measurement may be compared, in this 
sense, with the LEP measurements of ^4^^ (q=b,c). The latter depend, via Eqn(3.1), 
on Ai (defined by the values of the Z chargcd-lcpton couplings) and Aq (defined by the 
values of the Z heavy-quark couplings), Since Ai is independently measured, Aq can then 
be extracted from the measured value of A^/^. Similarly the NuTeV result depends, in 
a correlated way, on the values of sin 

Qon-sheii (equivalent to mw) and po (equivalent to 
the ZuV coupling). Since niw (and hence sin^ ^^"■''*^") is precisely determined at LEP 
and FERMILAB the correlated value of po, and so also s^,, can be extracted in a similar 
fashion to Aq from A'p^. 

The experimental values of the observables in Table 14 are compared with the SM 
prediction for mt = 178 GeV, mn = 120 GeV and a(mz) = 0.007755. This choice 
of mn (just above the experimental lower limit) is near to the maxima of the Ai{lept) 
and Ai{all) curves of logCL plotted in Figs. 2 and 3. Note that the model-independent 
analysis corresponds to only the Ai{lept) curves in Figs. 2-6. Only by making the stronger 
assumption of the SM values of the Z couplings to quarks, is it possible to derive Ai{had) 
and Ai(all). 

The agreement with the SM predictions shown in Table 14 is not completely satisfac- 
tory. The largest deviations are for Si, ( -0.84% and 3.7cr) and /l^ (-3.5% and 2.5cr). The 
positive 1.05(7 and 0.94(7 deviations of both Ai and mw respectively reflect the fact that 
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0.0774 
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-0.4209 




1996 
data 


9b 


0.1098(101) 


3.2 


[30] 


9k 


-0.4155(30) 


1.8 


1998 
data 


9^ 


0.1050(90) 


3.1 


[39] 


9k 


-0.4159(24) 


2.1 


2003 
data 


9b 


0.0951(63) 


2.8 


[4] 


9k 


-0.4182(16) 


1.7 



Table 16: History of measurements of and gk- The SM predictions correspond to rut 
= 178 GeV, ruH = 120 GeV and a{mz) = 0.007755. 



X 


1996 


2003 


[X(2003) 


- X(1996)]/ax(2003) 




0.0979(23) 


0.0997(16) 




1.1 


^b(LEP) 


0.3676(24) 


0.3663(13) 




-1.0 


^(SLC) 


0.863(49) 


0.925(20) 




3.1 


Ai{LEF) 


0.1466(33) 


0.1482(26) 




0.62 


AKSLC) 


0.1543(37) 


0.1513(21) 




-1.4 


A(LEP+SLC) 


0.1501(24) 


0.1501(16) 




0.0 



Table 17: Obscrvables in the 1996 and 2003 data sets contributing to measurements of 
the b-quark effective couphng constants. 
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the central values of mu preferred by these observables 

ruH = 72.51^^^ fit of Ai{lept) only 

niH — GS.Slgyl fit of mw only 
are incompatible with the direct lower hmit of mn — 114.4 GeV. 

The effective vector and axial vector couplings of charged leptons, neutrinos, c quarks 
and b quarks, that may be directly derived from the observables and Af {f = l,i>, Q) 
are presented in Table 15, in comparison with SM predictions. The agrement with the SM 
is satisfactory for the charged leptons and c quarks, but the neutrino couplings show a 
3.7(7, Vb a 2.5a and Ob a 2.9a deviation. As previously pointed out [16, 17] the apparently 
anomalous behaviour of the b-quark couplings is essentially found in the right-handed 
effective coupling, rather than the left-handed one, where: 

VS = ^ (7.4) 

t = ^ (7.5) 

It is interesting to consider the history of this apparent anomaly, which is illustrated 
in Table 16. The most significant deviation (42% and 3.2a) was seen in the 1996 data 
set [30]. In the current, essentially final, data set the size of the effect is reduced to 
32% but still has a significance of 2.8a. The left-handed coupling is now slightly more 
consistent with the SM (1.7o" deviation) as compared to 1998 {2.1a deviation) and 1996 
(1.8(7 deviation). The experimental error on g^ is reduced by ~ 40% in the current data 
set as compared to that of 1996. The sources of these changes are made clear by the 
entries of Table 17 in which are presented the model- independent observables used to 
calculate the b-quark couplings, as derived from the 1996 and 2003 data sets. Also shown 
are the shifts of the observables in units of the 2003 experimental errors. It can be seen 
that the most important change occurs in the direct measurement of Ai, from SLC The 
value of A;(LEP-|-SLC) used to extract Ab(LEP) from A^J^ is the same for the two data 
sets. 

To date, only a few authors [40, 41, 42] have proposed new physics interpretations of 
the measured b-quark couphngs. Other authors [13, 43] have argued that there is unlikely 
to be a new physics interpretation of the observed anomaly. The present writer finds all 
the reasons given for this conclusion to be either simply wrong, or unconvincing. It was 
argued that: 'the sensitivity of A^J^ to Ah is small, because Ai is small'. In fact, the size 
of Ai is irrelevant. What are important are the relative precisions with which it and A^J^ 
are known. Since the errors on Ai and A'p^ are uncorrelated it follows from Eqn(3.1) 
that: 

a(A,) 



\ 




(7.6) 



where a{X) is the experimental error on X. As shown in Tables 1 and 6, Ai and A^/^ have 
relative uncertainies of 1.2% and 1.6% respectively. Since the observed deviation from the 
SM is 3.6%, the contribution to the uncertainity on Aj, due to that on Ai is essentially 



^^These fitted values of mn are for a(mz) = 0.007755 and = 178 GeV 
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negligible in comparison with the the observed deviation from the SM prediction. It is 
further argued in Reference [13] that a new physics interpretation is disfavoured because 
a significant deviation is seen only in from LEP and not in Rb (equivalent to Sb) or 
Aii{SLC). But, in the case of a deviation from the SM only in the right-handed coupling, 
no significant change is expected in Rb^"^. Also the measured value of Ab{SLC) lies 1.5a 
above Ab(LEP) and O.Scr below the SM prediction. Furthermore XAb,wA/d-o.f. = 2.32/1, 
CL=0.13. There is therefore no strong evidence for any incosistency between the mea- 
sured values of ^^(SLC) and Af,(LEP). Both values are included in the weighted average 
that differs by 2.5a from the SM prediction. In Reference [13] it is stated that: 'One 
concludes that most probably the observed discrepancy is due to a large statistical fluc- 
tuation and/or an experimental problem'. The discussion in Section 3 above reaches just 
the opposite conclusion. In fact the statement just quoted tacitly implies (without any 
justification) that all the LEP experiments have seriously underestimated the systematic 
arrors of their A^p^ measurements; indeed, as discussed in Section 3 above, by at least an 
order of magnitude. This may be possible, but hardly seems likely. The total estimated 
experimental error on the LEP value of is is largely statistical. In contrast the Rb 
measurement has a statistical error of 0.20% to be compared with a systematic one of 
0.22%, so that the estimated systematic contribution is much more important than that 
for A^p^. Indeed from general experimental considerations, asymmetries such as ^^^^ 
are expected to have smaller systematic uncertainities than quantities such as Rb where 
absolute experimental detection efficencies play a role. Contrary to what is implied in 
Reference [13] then, there is no objective reason to suppose that the ^4^^ measurement 
should be less reliable than the Rb one. As discussed in Section 3, the hypothesis that the 
Ab deviation is a purely statistical effect has a CL of ^ 10~^ and so, though not completely 
excluded, is very unlikely. It remains true however that a correlated systematic error of 
unknown origin in the LEP A^p^ measurements is expected to produce an anomaly pre- 
dominantly in the right-handed effective couphng. Assuming SM values for the couphngs 
and that the -3.5% discrepancy in Ab is of systematic origin, the derived values of the cou- 
plings: g§ = 0.09486 and g{; = -0.4187 are in good agreement with the measured values 
in Table 16. Another argument [17] in favour of an unknown systematic origin for the Ab 
discrepancy is to note the good agreement of the measured value of Sb with the SM predic- 
tion shown in Table 14. This agreement requires the presence of of large, mf dependent, 
quantum corrections originating in the strong breaking of quark flavour symmetry in the 
third generation of SM fermions. Neglect of these corrections gives a prediction of 0.3707 
for Sb, differing from the measured value by 3.7cr. In the case of a new physics explanation 
of the Ab discrepancy, the appearence of the expected quantum corrections from the SM 
for Sb must be regarded as fortuitous. Thus, although there are no objective experimental 
reasons to doubt the correctness of the Ab measurement, a systematic effect of unknown 
origin cannot be excluded. The good agreement of Sb with the SM prediction and the 
large observed deviation in the right-handed coupling are consistent with this hypothesis. 
A purely statistical fluctuation is very unlikely. The effect could also be explained by new 
physics. There are no good reasons for the statement in Reference [13] that: 'It is well 
known that this (Ab) discrepancy is not likely to be explained by some new physics effect 
in the bbZ vertex'. In fact all the explanations mentioned above (including new physics) 
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r/iad(thy) definition 


Thad [GeV] 


[r/jad(thy)-r/jad(expt)]/cT(expt) 


SM 


1.7439 


-0.4 


b —>■ b(meas) 


1.7451 


0.35 


b,d b(meas) 


1.7414 


-1.50 


b,d,s b(meas) 


1.7377 


-3.4 




1.7444(20) 





Table 18: Constraints on the Z couplings to d-type quarks from the LEP average mea- 
surement of r had- 

remain open possibilities. Only better experimental data can decide between them. 

The next question that obviously arises is whether there is any evidence that the 
couplings of non-b quarks may also deviate from the SM predictions. Direct measurements 
of the light quark couplings have been performed by the DELPHI [44] and OPAL [45] 
Collaborations. For example, OPAL found: 

Tji _ _n 44+013 _ n 1 Q+0.15 

9d,s — > 9d,s — ^■^'^-0.17 

in good agreement with the SM predictions of -0.424 and 0.077 respectively. However, the 
very large uncertainties on the measurements preclude obtaining any useful information 
concerning deviations at the few % level, such as that observed for the parameter Af,. In 
fact the observed deviations of and Ac from the SM predictions by factors of 1.036 
and 1.023 are of comparable size. Since the relative error on Ac is two times larger than 
that on Ab, only for the latter is a possibly significant deviation from the SM observed. 
The similar qualitative behaviour of Ah and Ac with respect to the SM prediction can be 
seen in Fig. 15.1 of Reference [4]. The LEP measurement of can be used to extract 
an average value of Ag (averaged over all quark flavours) that is 1.07(7) times the SM 
prediction for this quantity. Thus, as previously mentioned, the different quark charge 
asymmetry measurements do not exclude deviations of Ag {q — u, d, s, c) as large as that 
observed for Ah- The present data are not, however, sufficiently precise to give any positive 
evidence for such an effect. 

As pointed out in Reference [46], much stronger constraints on the non-b quark cou- 
phngs are provided by the LEP average value of the hadronic width, Thad, of the Z 
boson [4]: 

Thad = 1.7444(20) GcV 

This value is in excellent agreement with the SM prediction^"'' of 1.7439 GeV (0.03 %, 
0.25(7 deviation). One may also note in Table 14 the almost perfect agreement of the s'^^ 
measurement with the SM prediction. The small relative uncertainty of 0.11% on Thad 
allows significant constraints to be placed on different hypotheses concerning the size of 
the Zqq couplings. Using the relation: 

/o/^ 3 u,d.s,c,b 

Tnad = ^^^^ E SgC^^^'C?^'' (7.7) 

"Since Sb '^[{'ShY + (ffb)^] and, with the SM values of the coupHngs: {g^f = 0.005, {g^f = 0.18, 
a 100% deviation ofTj^ from the SM prediction changes Sb (or Rb) by only 8%. 

^^Unless otherwise stated, all SM predictions are for nit = 178 GeV, niH = 120 GeV and a{mz) = 
0.007755. 
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the SM predictions for the d and s quarks may be replaced by the central values of the 
measured b-quark couphngs from Table 14. The results given by replacing, in Eqn(7.7), 

the SM predictions for (i) b quarks, (ii) b and d quarks and (iii) b, d and s quarks by the 
measured b-quark couplings from Table 14 are presented in Table 18, in comparison with 
the measured value of V^ad- It can be seen that, although the prediction is little changed 
for case (i), case (iii) is excluded by the measured value of Thad at the 3.4(T level. 

Since (see Table 14) the measured value of Sc agrees well with the SM prediction, the 
measured value of Vhad will provide no useful constraints if the procedure used in Table 
18 is repeated fou u-type quarks. The effective coupling constants and also agree 
well with the SM predictions. 

The experimental situation concerning measurements of right-handed and left-handed 
Z-fermion pair couplings and the W boson mass is summarised in Tables 19 and 20. In 
Table 19 the couplings of charged leptons, c quarks, b quarks and neutrinos are compared 
with SM predictions. Similar comparisons are made in Table 20, varying the values of 
nit and iTLu in the SM predictions. In this case, for clarity, only deviations from the SM 
predictions are tabulated. Tables 19 and 20 contain, in concise form, essentially all the 
precision information on the SM derived from the experimental programmes of LEP and 
SLD as well as the main contributions of FERMILAB to the same subject (essentially 
measurements of m^, mw and the ZvV coupling) during the same period. 

Looking at this comparison, it is difficult to conclude that the level of agreement with 
the SM is good. For mu = 120 GeV, i.e. around the maximum value of CL, as in Table 
19, show deviations of 2.8 and -3.7 standard deviations respectively. All of 

'gf', and niw show negative deviations around the one standard deviation level as a 
consequence of the low values of niH (inconsistent with the experimental direct lower 
limit) favoured by the measured values of these quantities. For mn = 300 GeV (see Table 
20) five out of the eight EW parameters show deviations around three standard deviations. 
Those associated with and are almost independent of uih and vary only weakly with 
mt- Another feature is that increasing (decreasing) the value of rrit improves (worsens) 
the agreement for and (itlw) for all values of mH- In any case, the SM still fares 
badly for itlh = 300 GeV and above. This is aleady apparent in Figs. 2 and 3, and will 
also be evident in the combined confidence level curves based on all precision data to be 
discussed below. 

Perhaps the most important aspect of the SM that has been tested by recent precision 
measurements is the renormalisability of the theory. This enables quantum loop correc- 
tions involving fermions, weak bosons and the Higgs boson to be calculated and compared 
to experiment, just as precise measurements of similar effects in QED, such as the Lamb 
shift of Hydrogen and the anomalous magnetic moments of the electron and muon, were 
important to establish the essential correctness of the theory for the description of such 
quantities. The model-independent observables shown in Table 14 can be used to isolate 
the effect of quantum corrections in the coupling of charged leptons, c quarks, b quarks 
and neutrinos to the Z boson. For this only the values of Ai{lept), Ac and Ai, are required 
for charged leptons and heavy quarks, and that of s^, for neutrinos. The A parameters are 
all simple mappings (see Eqn(3.2)) of the ratio: r — v/a of the vector and axial- vector 
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coupling constants. At tree-level in the SM the following relations hold: 



ri = 


1 - 




rc = 


1 - 




n = 


1 - 


3^^ 




1 






2 





r^' = 

Quant — 

= 

Quant 
Quant — 


(slY - 




l-ri 






3(1 -fe) 






3(1 -n) 






4«^ 


A^^ = 
Quant 




)/(i) = 


= 2^b - 1 



(7.8) 
(7.9) 

(7.10) 

(7.11) 

where Sy^ = 9'^~^^'^^^ as defined in Eqn(6.2). In the presence of quantum corrections 
Tf Tf, and {s'^Y (/ = l,c,b) in Eqns(7.8)-(7.10) and in Eqn(7.11). 

Thus parameters, ^Quanti ^^^^ measure directly the effect of quantum corrections can be 
introduced according to the definitions: 

(7.12) 

1 (7.13) 
1 (7.14) 
(7.15) 

In the absence of quantum corrections all the Sq^^^^ parameters vanish. The experi- 
mental values of these parameters derived from the entries of Table 14, as well as the 
corresponding SM predictions, are presented in Table 21. 

It can be seen from this table that the expected size of the corrections in the SM 
is ~ 4 — 5% for I, c and b, and an order of magnitude lower for u. By far the most 
significant measurement of quantum corrections (44 standard deviations from zero!) is 
that of Squanf Good agreement with the SM prediction {1.2a deviation) is found, at least 
for ruH = 120 GeV, as used in the SM predictions in Table 21. As discussed previously, in 
connection with and (see Table 20), the agreement worsens considerably for higher 
values of niH- For example, for niH = 300 GeV, (5q„„„((SM) = 0.04385, a deviation of 3.5a 
from the measured value. For c quarks the expected quantum correction is only a 1.7a 
effect. Good agreement with the SM prediction is found in this case. For b quarks the 
expected quantum correction is unmeasurable with present data (expected effect O.Gia 
from zero) whereas a S.Aa effect is actually observed, differing by 2.7a from the SM 
prediction. This is indicative, as discussed previously, of either a large and unknown 
correlated systematic effect (perhaps in combination with a statistical fluctuation) in the 
LEP A^p^ measurements, or of new physics at tree-level, or, indeed some combination 
of the two. A purely statistical fluctuation is very unlikely. The situation is exactly the 
opposite for the neutrino couplings. The SM predicts a large (4.2(j) quantum correction, 
while in the data only a 0.5a effect is seen. The discrepancy with the SM amounts to -3.7a 
in this case. Rather than a tree-level effect giving a large apparent quantum correction, 
as for the b quarks, it seems that the SM quantum corrections are effectively 'turned off' 
in the case of the the Zi/V couplings! Theoretical interpretations of this apparent coupling 
suppression have been made in References [47, 48] . 
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SM parameter 


Expt value 


SM 


(Exp-SM)/a 




0.23171(25) 


0.23202 


-1.2 


9i 


-0.26954(23) 


-0.26935 


-0.83 


9c 


-0.1585(48) 


-0.1547 


-0.79 


9c 


0.3460(36) 


0.3468 


-0.22 




0.0951(63) 


0.0774 


2.8 




-0.4182(16) 


-0.4209 


1.7 


t 


0.5003(6) 


0.50251 


-3.7 


mwlGeV] 


80.426(34) 


80.394 


0.94 



Table 19: Measured values of precision electroweak parameters compared to SM predic- 
tions for mt = 178 GeV, tuh = 120 GeV and a{mz) = 0.007755. 





rriH = 120 GeV 


rriH = 200 GeV 


mn = 300 GeV 


mt [GeV] 


173.8 178.0 182.2 


173.8 178.0 182.2 


173.8 178.0 182.2 


t 


-1.6 -1.2 -0.82 
-1.6 -0.83 0.03 


-2.4 -2.1 -1.7 
-3.0 -2.1 -1.3 


-3.1 -2.8 -2.4 
-4.0 -3.2 -2.4 


9c 

tc 


-0.77 -0.79 -0.80 
-0.19 -0.22 -0.27 


-0.75 -0.76 -0.77 
-0.11 -0.16 -0.2 


-0.72 -0.73 -0.75 
-0.06 -0.10 -0.15 


9^ 
9t 


2.81 2.82 2.82 
1.71 1.70 1.64 


2.80 2.80 2.81 
1.60 1.57 1.54 


2.79 2.79 2.80 
1.52 1.49 1.46 




-3.5 -3.7 -3.9 


-3.3 -3.5 -3.7 


-3.2 -3.4 -3.6 


mw 


1.7 0.94 0.18 


2.6 1.9 1.1 


3.4 2.7 1.9 



Table 20: Values of deviations (Expt-SM)/a for precision electroweak parameters. SM 
predictions with a{mz) = 0.007755. 

The main conclusion to be drawn from the results presented in Table 21 is that only 
for the charged leptons and c quarks is the present data reasonably consistent with the 
SM. This implies that since the c-quark couplings are almost completely insensitive to 
the values of rrit and itih (see Table 5) information on these parameters, via quantum 
corrections, can only be reliably obtained from the charged lepton couplings. This further 
implies (see Tables 11 and 12) that the maximum values of CL of 23% or 1.7% are 
obtained at mn = 115 GeV when the NuTeV mw measurement is excluded or included 
respectively. Much lower confidence levels of 0.054 and 0.0011 are obtained, under the 
same conditions, for mu = 220 GeV, the presently quoted 95% CL upper limit on uih 
from the EWWG [4]. In fact, even lower CLs will be found for a global analysis based on 
all model-independent observables shown in Table 14, which will now be discussed. 

In order to test the overall level of agreement of the precision data with the SM, 
estimators are calculated from the observables shown in Table 14. The corresponding 
confidence level CL(x|j^^) is then combined with that, CL^+b, of the direct Higgs search 
using Eqn(2.8) to give tuh confidence level curves, CL. Since all equivalent and statisi- 
cally compatible measurements are combined in order to extract the model-independent 
observables, there is, in this case, no contribution of the type Xx wa ^^^^ i^^o account 
the degree of statistical compatiblity of different measurements of the same quantity X. 
In view of the possibility of new physics or poorly understood systematic effects for the 
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[E.rpf - SM]/af 


/ 


0.04064(92) 


44.2 


0.04118 


44.8 


-0.59 


c 


0.060(25) 


2.4 


0.0418 


1.67 


0.76 


b 


0.249(74) 


3.4 


0.049 


0.65 


2.7 


V 


0.0006(12) 


0.5 


0.00502 


4.2 


-3.7 



Table 21: Quantum correction parameters for different fermion flavours f. SM predictions 
for mt = 178 GeV, mn = 120 GeV and a{mz) = 0.007755. 





NuTeV out 


NuTeV Si, meas. 


NuTeV m,w meas. 


mH (GeV) 








111 


4.9 X 10-^ 


4.9 X 10-« 


3.2 X 10-^ 


113 


2.0 X 10"^ 


2.1 X 10-5 


1.4 X 10-5 


115 


0.031 


3.8 X 10-3 


2.6 X 10-3 


140 


0.051 


5.1 X 10-3 


4.2 X 10-3 


180 


0.053 


5.9 X 10-3 


6.4 X 10-3 


220 


0.041 


5.6 X 10-3 


6.4 X 10-3 


260 


0.026 


2.9 X 10-3 


5.0 X 10-3 


300 


0.015 


1.6 X 10-3 


3.4 X 10-3 



Table 22: Combined confidence levels CL for consistency with the SM as a function of 
mH- AH data, as in Fig.6. 





NuTeV out 


NuTeV Sj^ meas. 


NuTeV mw meas. 


mu (GcV) 








111 


1.3 X 10-6 


7.4 X 10-^ 


4.7 X 10-s 


113 


4.5 X 10-^ 


2.8 X 10-5 


1.9 X 10-5 


115 


0.061 


4.7 X 10-3 


3.3 X 10-3 


140 


0.042 


2.2 X 10-3 


2.0 X 10-3 


180 


0.010 


5.3 X 10-^ 


7.4 X 10-^ 


220 


2.0 X 10-3 


1.0 X 10-^ 


2.1 X 10-"^ 


260 


3.5 X 10-^ 


1.9 X 10-5 


5.3 X 10-5 


300 


6.0 X 10-5 


3.4 X 10-6 


1.2 X 10-5 



Table 23: Combined confidence levels CL for consistency with the SM as a function of 
mH ■ Lepton data only, as in Fig. 7. 
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Figure 6: Combined mn confidence levels. CL^+f, is combined with CLIxsm) using 
Eqn(2.8). Xsm calculated from all model-independent variables: Ai{lept), si, Ac, A/,, 
Sb; SnbJ ^i^j and mw (NuTeV). 

Zbb couplings, and of the overall status, as well as the different possible interpretations 
of the NuTeV experiment, the above procedure is repeated for different selections and 
interpretations of the data: the results are presented in Tables 22 and 23 and Figs. 6 and 
7. 

In Table 22 and Fig. 6, all the observables in Table 14, except Sc, are included in the 
estimator. In the case that the NuTeV result is excluded, is assigned the LEP-only 
value of Table 5 and Eqn(6.7). The 'NuTeV meas.' curves use the Sj, value quoted 
in Table 14 and Eqn(6.9) (LEP+ NuTeV average). For the interpretation of the NuTeV 
result as an myy measurement, ?T74y(NuTcV) is included in Su is set to the LEP- 

only value. The number of degrees of freedom of the 'NuTeV out', 'NuTeV Si, meas.' and 
'NuTEV mw meas.' estimators are 10, 10 and 11 respectively. 

In the case of new physics, or unknown systematic effects, in the A^ measurement, it 
is clearly of interest to test only the level of agreement of the leptonic sector with the SM. 
For this, combined mn confidence level curves are also derived using only the leptonic 
observables: Ai{lept), si and as well as mw and, possibly mvF(NuTeV). The 'NuTeV 
out', 'NuTeV Su meas.' and 'NuTeV mw meas.' cases arc considered as previously, giving 
now 4, 4 and 5 degrees of freedom, respectively, for Xsm- The results for the 'Lepton data 
only' case are presented in Table 23 and Fig. 7. 

It can be seen from Table 22 and Fig. 6 that, when all data is included, the maximum 
of CL occurs around mn — 180 GeV, and has a value ~ 0.05 in the case that the 
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Figure 7: Combined itlh confidence levels. CLg+i i^ combined with CLIxsm) using 
Eqn(2.8). Xsm calculated from tie leptonic model-independent variables: Ai[lept), si 
'By, as well as mw and mw(NuTeV). 

NuTeV experiment is excluded and ~ 0.0065 when it is included, independently of the 
interpretation (si, or mw measurement) of the experiment. For mn = 300 GeV, the CLs 
are much lower, being 0.015, 0.0016 and 0.0034 for the 'NuTeV out','NuTeV meas.' 
and 'NuTeV mw meas.' cases respectively. In the case that only the leptonic Z-decay 
data is considered, the maximum value of CL occurs just above the direct lower limit 
with the values ~ 0.06 (NuTeV out) and ~ 0.004 (NuTeV included). For mn = 300 GeV, 
very low values of CL are found: 6.0 x 10~^, 3.4 x 10~^ and 1.2 x 10~^ for the three cases 
considered previously. 

8 Comparison with the EWWG and EWPDG Global 
Fits 

The conclusions presented above concerning the global level of agreement with the 
SM, and the possible value of mn may seem somewhat at variance with those drawn from 
the global fits of the EWWG [4] and EWPDG [5]. Indeed, the quoted confidence levels 
of the latter are much higher. The EWWG quotes confidence levels of 4.5% for 'all data' 
and 28% when NuTeV is excluded. These fits obtain (see Table 2) mH = QQ-ll GeV and 
mH — Qllle GeV respectively. These central fit values are completely inconsistent with 
the experimental direct lower limit on mn- Already for mn = 111 GeV, CLs+b — 10~^. 
In addition the fitted values of mn are strongly biased towards higher values, as discussed 



Lepton data only 




36 



above, due to the inclusion of the 'anomalous' quark asymmetry data in the fits. The 
latest global EWPDG fit, which did not include the NuTeV datum, found mn ^ QStH 
GcV and x^/rf.o./. = 47.3/38, CL= 14%. Some reasons for the higher confidence levels 
found for the global EWWG and EWPDG fits have been put forward in Section 2 above, 
so it is interesting, in the light of that discussion, to examine in detail the contributions 
of different types of observables to the of these fits. 

The 20 observables inclused in the EWWG fit may be classified as follows: 

Measured: A|f2^(m|), mz, rrit. 

m^^-sensitive: Ai{P^), sin^ Q'^Jj{Q''/^), Ai{SLD),Afs, A%,mw,sm^ Qw{l^^f)■ 

OtheJ:s■. alTz, Rl RlRlAhAc ^w, Qw{Cs). 

The relative contributions to the total of the different types of observable are 
summarised in the first three rows of Table 24. The quantities I] A(x')^ are calculated 
from the squares of the 'pulls' given in the last column of Table 16.1 of Reference [4]. Also 
shown in Table 18 are the 'pseudo-confidence levels', CL', for each A(x')^ assuming that 
the effective number of degrees of freedom, d.o.f.', is equal to the number of observables 
of each type. Since there are five fitted parameters: mz, Mt, log(mij), a{mz) and as{mz), 
the true total number of degrees of freedom is 15 rather than 20, so that the quoted values 
of CL' should be considered as upper limits on the true CL to be associated with each 
type of observable. This procedure neglects correlations between the fitted parameters, 
but is adequate to show the very different contributions of the different observable types 
to the overall x^. The 'Measured' observables, i.e. those that are identical to fitted 
parameters, are seen to provide an essentially vanishing contribution to x^, since the 
fitted parameters are completely determined by the corresponding observables. Thus 
three extra degrees of freedom are obtained, free of charge, which considerably improve 
the CL of the SM comparison^^. The Measured parameters should be fixed in the fit, 
not treated as observables to be fitted. The effective confidence level, CL' of the 8 'mjy- 
sensitive' observables is a factor of 21 smaller than the CL' of all observables (see Table 
24) and a factor 6.4 smaller than the EWWG global fit CL. The value of CL' for the 
m/j-sensitive observables is similar to the maximum value 0.0064 of CL for the analysis 
of all model-independent observables in the last column of Table 22 and Fig.6. The 
'Other' observables give a rather low contribution to the total (x')^ : as compared to the 
expectation from the number of degrees of freedom, which give a further improvement to 
the CL of the global fit beyond that expected from the inclusion of observables that are 
only weakly sensitive to rrit and niH- Excluding the NuTeV datum sin^ 0vk(^A/') from the 
group of m//-sensitive observables gives A(x')^ = 12.6, d.o.f = 7, CL'= 0.08 which 
is similar to the maximum value of CL of 0.053 for the corresponding analysis of model- 
independent observables in Fig 6 and the first column of Table 22^^ Thus the difference 

^^For example, calculating CL' by summing the values of A(x')^ corresponding to each observable 
type gives (x')"^ /d.o.f .' = 26.64/20, CL'= 0.15 when the Measured observables are included (see Table 
24) and {x'f /d.o.f.' = 26.60/17, CL'= 0.064 when they are excluded. 

i^Actually the EWWG and EWPDG fits use the old mt value 174.3(5.2) GeV rather thanJJ8.0(4.3) 
as used in Tables 22 and 23. This, however, has only a minor effect on the maximum values of CL. Using 
the old value of TOj, maximum values of CL of 0.0080 [0.057] are found when the NuTeV experiment is 
included [excluded], to be compared with the corresponding numbers 0.0064 [0.053] quoted above using 
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between the confidence levels found in the present paper and those of the global EWWG 
fits can be largely understood as a consequence of the extra degrees of freedom associated 
with the 'Measured' and 'Other' observables which are only weakly sensitive to the crucial 
unknown parameter, ttih, of the SM. 

Classifying the 42 observables used in the latest global EWPDG fit [5] in the same 
fashion as done above for the EWWG fit gives: 

Measured: mz, rrit. 

m/z-sensitive: Ap^, Ap'^, Ap^, Ap^, Ap^, Ap^,sm'^ Q^fj:{Q'^^),Ae{l),Ae{2),Ae{S), 
Af,,Ar{l),Ar{2), miy(FERMILAB), miy(LEP) 

Others: Fz, T^ad, '^inv, ^l+l-, 0"°, -Rg, -Rfe,-Rc,^6,^c,^s, R~ iH'^-, -R'^(l), R'^{2), 

g'(.'{l),g'^%2), g\-{l),g\-{2), Qw{Cs),Qw{Tl),T{b ^ sj)/T{c ^ eu), {g^-2-a/n)/2 

The difference with respect to the EWWG fit is that a wider range of observables are 
included as well as several different measurements of the same quantity, indicated, for 

example, as Ae{l),Ae{2)..., and that charged lepton universality has not been used to 
reduce the number of observables. It could be argued that some of the 'Other' observables 
such as Fz, Thad and Fi„^ are actually quite 'mj^-sensitive', but the choice of the latter 
type of observable has been restricted, for purposes of comparison, to correspond as closely 
as possible to that made above for the EWWG fit. The values of 5] A(x')^; d.o.fJ and 
CL' for the three classes of observables in the EWPDG fit are presented in the fifth, 
sixth and seventh rows of Table 24. There are only two Measured observables, rrit and 
iTLz-, since the fitted parameters are: mf, m^, itlh and a{mz). Unlike for the EWWG 
fit as{mz) is treated as a fixed rather than a fitted parameter. The high value of CL' 
for the Measured observables shows that, as in the case of the EWWG fit, it is more 
appropriate to treat rrit and mz as fixed parameters in the fit. As for the EWWG fit, the 
value of CL' for the m/^-sensitive observables is much less than the global fit confidence 
level. The value of CL' for the 'Other' observables gives no indication for a possible over- 
estimation of systematic errors as in the EWWG fit. However, the observables Fj„„ and 
(g^ — 2 — a/7r)/2 not included in the EWWG fit show quite large deviations from the 
SM prediction. Removing these observables from the 'Others' set gives: I] A(x')^ = 14.4, 
d.o.f.' = 23, CL'= 0.91, again indicating a possible overestimation of systematic errors 
for the remaining observables, leading to a confidence level for the global fit that gives an 
optimistic estimate of the level of agreement with the SM. The lower confidence level of 
the global EWPDG fit as compared to the EWWG one, in spite of the inclusion of many 
unaveraged observables, is also explained by the inclusion of these 'discrepant' observables 
in the former fit. Excluding them gives x^/ d.o.f. ~ 37.7/36, CL = 0.39 which is larger 
than that of the EWWG fit (0.28), as expected. 

Thus the lower confidence levels found in the global analysis of the present paper 
as compared to those quoted by the EWWG and EWPDG are fully explained in terms 
of the dilution of the hypothesis testing power of the latter fits due to the inclusion of 
unaveraged or insensitive observables in the estimator. 

the new rrit measurement. 
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Observable Type 




d.o.f.' 


CL' 




Measured 


0.04 


3 


0.998 


EWWG 


mif-sensitive 


21.0 


8 


0.0071 




Others 


5.6 


9 


0.78 




All 


26.6 


20 


0.15 




Measured 


0.05 


2 


0.975 


EWPDG 


TOff-sensitive 


25.0 


15 


0.050 




Others 


25.1 


25 


0.78 




All 


50.2 


42 


0.18 



Table 24: Contributions of different types of observables to the of the latest global 
EWWG [4] and EWPDG [5] fits. See text for definitions of the quantities shown. 

A final remark concerning the CLs quoted for the EWWG and EWPDG fits is that, as 
previously mentioned, these numbers correspond to central fitted values of mn that are 
grossly inconsistent with direct the experimental lower limit of 114.4 GeV. Replacing the 
fitted values by this limit will evidently yield lower CLs. Because of the relatively large 
uncertainy on the fitted value of mn, the expected reduction in the CL will be less than 
that due to correcting for the statistical dilution effects discussed above. It still remains 
however a meaningless exercise to quote CLs for fits with mn — 90 — 100 GeV, values 
that are excluded by the direct search result with a CL (see Fig.l) very much less than 
10-7. 

The latest EWWG report [4], quotes a 95 % confidence level upper limit on mn of 219 
GeV. This estimate is based on the famous 'blue band'^^ plot (Fig 16.5 of Reference [4]) 
which shows the quantity: Ax^ = x'^~Xmin ^'^^ global fit as a function of rriH- Choosing 
an appropriate fixed value of Ax^ to define a confidence limit for rriH is a valid procedure 
only in the case that the confidence level derived from the value of of the fit to the 
mi/-sensitive variables is sufficiently high. Inspection of Figs 2 and 3 and Tables 11 and 
12 shows that this is hardly the case for the actual electroweak data set. It was argued in 
the previous Section that most reliable estimate of itlh is that derived using only charged 
lepton asymmetry data and mw- In this case the the maximum value of CL is 0.017 when 
the NuTev mw measurement is included, 0.23, if it is excluded. Thus the A^^ estimator 
for niH is acceptable only for the case of a fit to Ai{lept) and mw excluding the NuTeV 
measurement. However, at the 95 % confidence level EWWG upper limit of mn = 220 
GeV, the corresponding values of CL are 1.1 x lO^'^ and 0.054 with the implication that 
if a Higgs boson existed with this mass, and including, as in the final EWWG fit, the 
NuTeV measurement, the SM would be excluded by the data with a confidence level only 
slightly larger than 10""^! This is not at all the message that one might naively draw 
from the corresponding EWWG numbers. A confidence level of 4.5% for a global fit gives 
the impression that the data is, in general, not too badly described by the SM. In fact 
as a brief inspection of Tables 19-23, which contain all relevant information, show, this 
is hardly the case, so that a meaningful estimate of mn cannot be derived from a Ax^ 
plot based on such a global fit. In the model-independent analysis of all precision data 
in the previous Section, the level of the discrepancy with the SM may be even larger 

So-called becuse of the coloured version shown in dozens of electroweak review talks over the last 
decade 
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than if only the mif-sensitive observables Ai and mw are considered. Referring to Figs 
6 and 7 and Tables 22 and 23 and taking the ' NuTeV s^, meas.' curves, mandatory in 
a model-independent analysis, gives values of CL of 0.0038 (0.0056) for rriH = 115 (220) 
GeV for the 'All data' case and 0.0047 (1.0 x 10^'*) for the most reliable 'Lepton only' 
data. The reader must judge for herself (or himself) whether, in these circumstances, the 
SM does, or does not, provide an adequate description of the current data. Only in the 
case that it does, the 'blue band' plot, derived from a global fit to this data, provides a 
meaningful upper limit on niH- 



9 Summary and Conclusions 



It has been demonstrated in this paper that the the statistical estimator, XdataSM 
universally employed by the EWWG and EWPDG to judge the level of overall agreement 
of precision electroweak data with the SM predictions typically yields a confidence level 
an order of magnitude higher than estimators chosen to test specifically this level of 
agreement rather than the internal consistency of different measurements. The reasons for 
this are discussed in Section 2. While, as shown in Section 4, the number of independent 
observables sensitive to the most important poorly known parameters rrit and mn is 
very small, large numbers of observables (20 for EWWG, 42 for EWPDG) are used to 
construct Xdata,SM- these circumstances the corresponding CL reflects more the internal 
consistency of measurements of different observables than the level of agreement of the 
essential 'refined' parameters with the SM. Further dilution of the hypothesis testing power 
of Xdata SM Tcsults from the inclusion of many observables only weakly sensitive to rrit and 
iJiH as well as fit parameters identical to measured quantities that give anomalously low 
contributions to the global x^. 

In Section 3 the internal consistency of different heavy quark asymmetry measurements 
is discussed in detail. It is found to be very good. There is a hint, from the data itself, that 
uncorrelated systematic errors may be overestimated. Thus the observed 2.5cr deviation 
of the LEP+SLD average value of Aj, may be an underestimation. Combination of the 
statistically independent ^-^d Run Tests for the LEP and SLD measurements of Aj, 
yields a CL of 2.1 x 10~^ for a pure statistical fiuctuation, which therefore seems unlikely. 
Attributing the Ab deviation to a correlated systematic error of unknown origin in the 
LEP A^p^ measurements would require an effect that is 1.8 times larger than the estimated 
QCD correction, which is the dominant source of correlated systematic uncertainty on 
this quantity, and 13 times larger than the estimated uncertainty on this correction. 
Even given the inevitable theoretical uncertainties associated with QCD effects this again 
seems unlikely. Thus there is no sound experimental reason to doubt, as suggested by some 
authors [13, 43], the validity of the b-quark asymmetry measurements and the possible 
evidence they provide for new physics beyond the SM. 

It is shown in Section 4 that essentially all information from quantum loop effects on 
the values of rrit and uih is provided by only three observables: Ai{lept), si and niw (see 
Table 5). For ttih there are only two strongly sensitive observables: Ai(lept) and mw It is 
demonstrated by fitting only the mij-sensitive observables and also complete sets of model- 
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independent observables that indeed the fitted value of uih is essentially determined by 
the former set only. In an approach similar to that followed in Reference [19] it is assumed, 
in Section 4, that the SM is valid so that different values of Ai may be extracted from 
purely leptonic data {Ai{lept)) and from quark forward/backward asymmetries {Ai{had)). 
As previously noted [19], for the case of the equivalent observable sin^B^J^^, these two 
estimates differ by 3. Oct showing that the SM interpretation of the data is also inconsistent 
at this level. This is a simple consequence of the observed anomalous behaviour of Ai, 
and the possibly similar behaviour of other Zqq couplings. Also, as previously noticed, 
much larger values of itlh are favoured by Ai{had) than by Ai(lept). This is because 
(see Table 5) in the SM only Ai, not A^, is sensitive to mn and the measured forward 
backward asmmetry A^p^^ is, as shown in Eqn(3.1), proportional to AiAi,. Whatever the 
explanation of the A}, anomaly, the most reliable estimate of mn that can be derived from 
the present data is therefore that provided by Aiilept). Inclusion of Ai{had), assuming 
the correctness of the SM, gives a calculable bias towards higher values of uih resulting 
from the Ai-Ai, correlation in the quark forward/backward asymmetry. 

In Section 5 a combined confidence level, CL, as a function of mn is derived from the 
directly measured [10] confidence level curve CLg+t and the CL curve derived from the 

estimator using the m//-sensitive observables Ai and niw (see Fig.l). As in Section 4 
the correctness of the SM is assumed and different CL curves are calculated for Ai(lept), 
Ai{had) and Ai{all), where the latter is the weighted average of the former two observables. 
The above mentioned inconsistency between Ai{lept) and Ai{had) is taken into account 
when calculating the CL curve for Ai{all). Consistent results for CL (see Tables 9 and 
10) are found using either Xdata,wA + Xwa,sm directly or combining the CLs of Xdata,wA 
and XwA,SM using Eqn(2.8). Only the standard interpretation of the NuTeV experiment, 
as a measurement of mw, is considered in Section 4, and the CL curves are calculated 
both including and excluding this datum. The results for CL are shown in Tables 9-12 
and Figs. 2 and 3. The inconsistency of the SM interpretation is evident on inspection 
of these figures. For ttih = 300 GeV values of CL differing by two orders of magnitude 
are obtained from Ai{lept) and Ai{had). Also studied in Section 5 is the dependence 
of the CL curves on the assumed values of rrit and a{mz) (Table 13 and Figs. 4 and 5). 
Variation of by plus or minus the experimental uncertainty changes the value of CL by 
more than three orders of magnitude for rriH — 300 GcV; somewhat smaller changes are 
given by a similar variation of a(mz). This demonstrates the importance of more precise 
measurements of these parameters in order to obtain well defined SM predictions. 

In Section 6 the alternative interpretation of the result of the NuTeV experiment as a 
measurement of the ZuV coupling, rather than mw, is considered. It is pointed out that 
the former interpretation, mandatory in a model-independent analysis where no a priori 
assumptions are made concerning the strengths of the Zff couplings, is highly favoured 
by an argument based on the internal consistency of LEP and NuTeV data. In this case 
the model-independent observable s^, derived from the LEP and NuTeV data differs from 
the SM prediction by 3.7(t and so is the largest single deviation observed from a SM 
prediction. The combined CL, taking into account both agreement with the SM and data 
consistency, is very similar for either interpretation of the NuTeV result. 

An analysis in terms of model-independent observables similar to those previously 
published for earlier electroweak data sets [16, 17, 18, 46] is presented in Section 7. Results 
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are presented in terms of the 'maximally uncorrelated' observables presented in Table 14. 
Vector and axial vector couplings of the Z to fermion pairs (Table 15), and the equivalent 
right-handed and left-handed couplings (Table 19) are also presented and compared with 
SM predictions. The history of, and the different possible physical interpretations of, the 
Af, anomaly are also discussed. Although there are no purely experimental reasons for 
doubting the correctness of the fully compatible LEP and SLD measurements of At, it 
is pointed out that the good agreement between the measured values of Sb and the SM 
prediction (which requires the presence of large nit dependent quantum corrections) must 
be fortuitous if the Ah anomaly is to be explained by new physics. This is an argument 
suggesting an unknown systematic bias as the cause of the effect. It is also pointed out 
that such a systematic bias would result in an anomaly predominantly in the right handed 
couphng (2.8 a effect observed) rather than in the left-handed one (1.7 a effect observed). 
This implies that only new, more precise, experiments can discriminate between new 
physics and unknown systematic bias as the cause of the observed anomaly in the right 
handed b-quark coupling. It is shown that the measured value of r^a,^ renders unlikely 
the possibhty that the couplings of the other d-type quarks differ from the SM prediction 
in the same way as observed for the b-quark couplings. The existing direct measurements 
of these couplings are insufficiently precise to provide any useful constraints. 

In Section 7 the observed and expected quantum loop corrections are also discussed 
(see Table 21). The most significant and the most precisely measured effect, 5 Quanta defined 
in Eqn(7.12) is a 4% effect measured with a relative precision of 2.2% (measured effect 44(7 
from zero). For a value of niH of 120 GeV, the agreement of Sq^^^^^^ with the SM prediction 
is quite satisfactory (1.2o" deviation) but, as shown in Table 20, for higher values of itlh 
the level of agreement of the charged lepton couplings that determine Sguant deteriorates 
rapidly. Sg^^^^ (Eqn(7.13) is also predicted to be ~ 0.04 in good agreeement with the 
measurement (0.7cr deviation) but here the fractional precision of the measurement is 
only 42%. is predicted to be ^ 0.05 with an expected relative accuracy of 151%, 

so that no significant measurement is to be expected in this case. In fact the measured 
value of Sq^^^^ is much larger, 0.249(74), and so requires new physics at the tree level 
(2.7(7 deviation from the SM) if the data is correct. In contrast the experimental value of 
^Quant (Eqn(7.15) is 0.00502 with an expected relative experimental uncertainty of 24% 
(4.2(7 deviation from zero) whereas the measured value is 0.0006(12) (0.5ct deviation from 
zero). In this case the expected quantum corrections are not observed, leading to a -3.7ct 
deviation from SM prediction. Some theoretical interpretations of this effect have already 
been proposed [47, 48]. 

Finally, in Section 7, combined confidence level curves CL are derived including in 
the estimator not only the mj^-sensitive observables, as in Section 5, but all, or cho- 
sen subsets of, other model-independent observables. In order to see the impact of the 
NuTeV experiment, CL curves are are presented excluding the NuTeV data or for the 
two alternative interpretations: s,y or mw measurements. The same set of CL curves is 
also obtained using only the leptonic observables: Ai{lept), si and s,y in addition to mw- 
The results are shown in Tables 6 and 7 and Tables 22 and 23. The curves where NuTeV 
is included lie about an order of magnitude below those where it is excluded. The CL 
curves for the two different interpretations of the NuTeV experiment are similar, except 
that those corresponding to an Si, measurement lie significantly lower for large values of 
ruH of > 300 GeV. 
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Figure 8: Combined mn confidence levels. See Fig.6 for the definition of 'All data'. 



In Section 8 a comparison is made between the CLs found previously in the present 
paper and those quoted for the latest published EWWG and EWPDG global fits. They 
arc shown to be quite consistent when the various dilution effects of the estimators 
used in the global fits are taken into account. The maximum value of CL using all data 
and choosing the measurement interpretation of the NuTeV experiment of 0.0059, at 
TTiH = 180 GeV, is about an order of magnitude lower than the CL of 4.5% quoted for 
the 'all data' EWWG fit. The CLs of the global EWWG and EWPDG fits correspond to 
central fitted values of itlh completely excluded by the direct experimental lower limit of 
114.4 GeV at 95% CL. Replacing the fitted values of niH by this lower limit, to give the 
largest possible CL consistent with all experimental data, will result in lower CLs than 
those quoted for the fits. 

Finally are shown, in Fig. 8, the author's personal choice of the three most pertinent 
CL curves among the 18 different ones previously presented in this paper. The 'Ai{lept) 
and mw only' curve (dotted) gives the most rehable estimate of mn- The value of CL 
of ~ 0.2 — 0.3 for values of uih just above the direct lower limit of 114.4 GeV is quite 
acceptable. However for ttih = 300 GeV, CL is < 10~^ implying that, if the SM describes 
correctly the charged lepton sector and the Higgs boson exists, it must be very light indeed: 
< 180 GeV if CL > 0.05. Higher values of tuh are favoured by the 'AH data NuTeV out' 
curve (dashed). This is mainly due to the Ai — A^ correlation following from Eqn(3,l) 
in the measurement, as discussed above. The maximum value of CL is ~ 0.05 at 
about iTiH = 140 GeV. Including the NuTeV measurement gives the 'All data NuTeV s^, 
meas.' curve (solid line) with a similar shape but lying roughly an order of magnitude 
lower. The maximum value of CL is ~ 0.006 at uih — 180 GeV. This accurately reflects 
the best possible level of agreement, with the SM prediction, of the entire electroweak 
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data set. It is an order of magnitude, or more, lower than the CLs quoted by the EWWG 
and EWPDG groups for their global fits to similar data sets. 
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